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I. INTRODUCTION 



In the last decade we have witnessed significant advances in the field of trapped ultracold gases [2] opening new av- 
enues for the investigation of low-dimensional physical systems which can be well approximated by integrable models. 
The paradigmatic example is the Bose gas with contact interaction [3J, also known as the Lieb-Liniger model, whose 
experimental realization [2HS] has spurred renewed interest in computing physical properties which are experimentally 
accessible. In particular, the correlation functions, which can be measured using interference |10H13) . analysis of par- 
ticle losses pj Q3] , photoassociation [T3] , Bragg and photoemission spectroscopy [H)H2"U] , density fluctuation statistics 
[HJ[22], time-of-flight correlation statistics [33J and scanning electron microscopy [21] are extremely important. 

Despite the integrability of the model the calculation of the correlation functions is an extremely challenging problem 
which remains unsolved to this day. Significant simplifications occur in the case of infinite repulsion when the system 
is equivalent to free fermions. In this case the correlators can be expressed as Fredholm or Toeplitz determinants 
and the asymptotic behavior can be extracted from the solution of an associated Riemann-Hilbert problem [25H33) . 
Similar results, albeit in a non-rigorous fashion, can be derived using the replica method [34] . 

The introduction of the algebraic Bethe ansatz (ABA) provided the necessary tools to tackle the harder problem 
of calculating the correlation functions of integrable models away from the free fermion point [551458] . At zero 
temperature, the members of the Lyon group (Kitanine, Kozlowski, Maillet, Slavnov and Terras), making use of the 
results obtained in [3"9"H4"3"] derived in [H] the asymptotic behavior of the static density correlators in the repulsive 
Lieb-Liniger model and similar results for the longitudinal correlation of the XXZ spin chain. The large-distance and 
long-time asymptotic analysis of the density-density and field-field correlators was performed in 45 47 . In all cases 
these exact results reproduce and generalize the predictions of the Tomonaga-Luttinger liquid(TLL)/Conformal Field 
Theory (CFT) approach [48 -50 . A method of determining the "non- universal" prefactors appearing in the TLL/CFT 
expansion was introduced in [51] [52] and, also, very recently, in [53] . 

The temperature dependent correlation functions of a system characterized by a Hamiltonian H arc defined as 

\°)T- £ e _ B / T ' W 

where O is a local operator, the sum is over all the eigenstates |f2) of the Hamiltonian, and E their respective 
energies. The summation appearing in ([!]) makes the calculation of temperature dependent correlation functions 
extremely difficult. However, in the case of the interacting Bose gas we can circumvent this problem in two ways. 
First, it can be shown, see Chap. I of [35] and references therein, that in the thermodynamic limit can be replaced 

by 

{<d)t ~ (W (2) 

where \Qt) is any of the eigenstates corresponding to thermal equilibrium. This allowed the authors of 54, 55] to 
employ a method similar with the zero temperature analysis performed in [44] to obtain the asymptotic expansion 
of the generating functional of density correlators. The second method utilizes the quantum transfer matrix (QTM) 
and the connection between the XXZ spin chain and the one-dimensional Bose gas. Introduced and developed in 
[561 157 ) . the QTM, in particular its spectrum, is an extremely important tool in the investigation of temperature 
dependent properties of lattice systems. The free energy of the system is related to the largest eigenvalue of the 
QTM [53 [55] and the correlation lengths of the Green's functions can be obtained as ratios of the largest and next- 
largest eigenvalues [58] [59] . At the same time the QTM is a fundamental ingredient in obtaining multiple integral 
representations for temperature dependent correlation functions 60 62 . Even though there is no QTM equivalent 
for continuous systems we can use the fact that the one-dimensional Bose gas can be obtained in a specific continuum 
limit of the XXZ spin-chain [S3] [M] ■ Performing this continuum limit in the non- linear integral equations (NLIEs) 
characterizing the eigenvalues of the XXZ spin-chain QTM we will obtain the spectrum of what we will call the 
"continuum" QTM, from which we can calculate the thermodynamics and the correlation lengths of the continuum 
system. It is precisely this method that we will use in this paper to obtain the asymptotic expansion of the temperature 
dependent density-density and field-field correlation functions in the one-dimensional Bose gas. We should mention 
that the same scaling limit was used in |65j to obtain the k-body local correlators, i.e., correlation functions of the 
type ((5'^(0))' c (\E , (0))' c )t, for k < 4 using, however, a totally different method than ours. For k < 3 local correlators 
were first calculated in [B1>HT0"] . Other important results concerning the correlation functions of the ID Bose gas can 
be found in HH 

The plan of the paper is as follows. In the next section we introduce the one-dimensional Bose gas and present the 
asymptotic expansions for the correlation functions which constitute the main results of this paper. In Sec. |III| we 
review the XXZ spin chain and introduce the continuum limit which allows for the derivation of the Bose gas results. 



3 



In Sec. IV we introduce the XXZ spin chain QTM and obtain NLIEs for the largest and next-largest eigenvalues from 
which the correlation lengths can be extracted. The validity of the asymptotic expansions is checked in Sec. [V] by 
comparison with the TLL/CFT predictions. Finally, the asymptotic behavior of the correlators in the Bose gas is 
obtained by taking the continuum limit in Sec. |VI| Some technical calculations are presented in several appendices. 



II. THE ONE-DIMENSIONAL BOSE GAS AND MAIN RESULT 



We consider a one-dimensional system of bosons interacting via a <5-function potential with periodic boundary 
conditions. The relevant Hamiltonian is 

H NLS = [ dx [d x ¥(x)d x ^(x) + c^(x)¥(x)^{x)^(x) - n&(x)*(x)] , (3) 
Jo 

where c > is the coupling constant, \x the chemical potential, I the length of the system and we have considered 
h = 2m = 1, with m the mass of the particles. In ^ ^ (x) and ty(x) are Bose fields satisfying the canonical 
commutation relations 

¥(x')} = S(x - x') , [¥(3), (a/)] = [&(x), &(x')\ = . 

The interacting one-dimensional Bose gas, also known as the Lieb-Liniger or the quantum Non-Linear Schrddinger 
(NLS) model, is solvable by Bethe ansatz [3J[3E1|SD]. In the case of n particles the energy spectrum is given by 

n 

^(W) =!>(**), e (fc) = fc 2 - M , (4) 

j'=l 

with the quasimomenta kj satisfying the following set of Bethe ansatz equations (BAEs) 

<° ii ill: ;=i,-.». w 

Iv i rug fcU- 

It is useful to present the logarithmic form of the BAEs (J5j) 

n 

Ikj + 0(kj — k s ) = 2irrrij , 

s=l 

where rrij are integers or half-integers and the scattering phase 9{k) is defined by 

0(k) = Hog ( t -^~] , lim 9(k) = ±tt . 

\ic~ k J k^±oo 

At zero temperature and fixed number of particles n the ground state is obtained when the (half) integers take the 
values rrij = j — (n + l)/2 , j = 1, • • • , n [3]. In the thermodynamic limit I, n — > oo , with their ratio finite D = n/l, 
the values of the momenta kj condense in the interval [— q, q] called the Fermi zone or Dirac sea and the following 
integral equation for the density of particles in momentum space can be derived: 

m - I j\K(k - W) *' = I , K(k - k>) = P(k - k') = {k J; 2 + c2 • (6) 

The Fermi momentum q can be obtained as a unique function of D, the density of particles, via D = n/l = p(k) dk . 

At finite temperature the thermodynamics of the model was calculated in [T] (for a rigorous derivation, see [§Tll92| ). 
The grandcanonical potential per length is given by 




with s(k), the dressed energy, satisfying the Yang- Yang equation 

e(k) = k 2 - ii - — f K{k- k') log f 1 + e - wt - k '^ T ) dk' . (8) 

27T J M V I 
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A. Main result 



The main result of this paper is the computation of the large-distance asymptotic behavior of the correlation 
functions in the ID Bose gas at finite temperature. Due to the fact that the derivation of the asymptotic expansions 
is quite involved we prefer to present these results in the beginning of the paper. The interested reader can find the 
details in the following sections. 

We will start with the static density-density correlation function, (j(x)j(0))T, with j(x) = w{x)^{x). Consider 
the following set of functions Uj(fc) satisfying the nonlinear integral equations 

u l (k) = k 2 - fi + iT^6(k-k+)-iT^T6(k~k-)- |- J K(k - k ')log (l + e -^ (k ' )/T ^ dk! . (9) 

The 2r parameters, {k^} T j=i appearing in Eq. |9| are located in the upper (lower) half of the complex plane 

and satisfy the constraint 

1 + e ^(/c±)/T = a 

For a given r, the previous equation has more than 2r solutions, the subscript i labels all the possible choices of 
solutions for all r = 1,2, • • • . Note that the NLIEs ^ are almost identical with the Yang- Yang equation for the 
dressed energy ^ with the exception of the additional driving terms. The large distance asymptotic expansion for 
the density-density correlation function has the form 

(j{x)j(0)) T = const + M e~« (d) ^i , x -> oo , (10) 

i 

where are distance independent amplitudes which cannot be obtained using our method and the correlation lengths 
are given by 



1 If, /l + e -«*W/ T 



, , - , log ( - ' - _,, w ^ ) dk - kf +iS^kJ ■ (11) 



with e(k) the dressed energy satisfying (18]). Comparison with the TLL/CFT expansion (100) and other exact results 



(Chap. XVII of [55]) allows the identification of the constant term with (i(0))y. The leading terms in the expansion 



(10) are obtained considering r = 1 in Eq. km with the parameters kf , satisfying \-\-e Ui ^ kl ^ T — 0, closest to the real 



axis. 

A few remarks are in order. Using a different method almost identical equations were obtained by Kozlowski, 
Maillet and Slavnov [54] [55] for the generating functional of density correlators, (e^Jo ^ x ^ dx )t, from which the 

density correlator can be obtained via {j{x)j{Q))T — \-§^i-§^i {e v ^> ^ x ^ dx ) T 1 The only difference between 

our equations and the ones derived in [541 155] is the presence of a renormalized chemical potential fj, — > ji + <pT in 
the r.h.s. of Eq. As we will show in Appendix [C] a slight modification of our method allows for the derivation of 
the asymptotic expansion for the generating functional. However, in order to not confuse the reader, we prefer here 
and in the following sections to focus on the density and field correlators (see below) because it allows for an almost 
similar treatment. 

In the case of the field- field correlation function (vf^ (x)\&(0))t we introduce the set of functions Vi(k) satisfying the 
NLIEs 

v t (k) = k 2 - n±inT + iTy^9{k - k+) - iTV(9(fc - fcr) - — / K(k - k') log f 1 + e -M k ')/ T ) dk 1 . (12) 

3=0 3=1 JR 

The functions Vi{k) depend on 2r + l parameters, {k^} T j=o ; ({^7}J=i) located in the upper (lower) half of the complex 
plane which satisfy the constraint 

l + e ^ ± )/r = . 



It should be remarked that the authors of |54l I55| noticed that their results which were derived using the asymptotic analysis of a 
generalized sine-kernel Fredholm determinant can be interpreted in the framework of the QTM which is the primary object of this paper 
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In Eq. ( 12 1 we will consider the plus sign in front of the ittT term when is in the first quadrant of the complex plane 



5Rfc+ > , > 0, and the minus sign when k$ is in the second quadrant of the complex plane Kfeg < , Qk^ > 0. 
As in the case of the functions Uj(fc) the subscript i labels all the possible choices of roots for all r = 0, 1, 2, • 
large distance asymptotic expansion of the field-field correlation function has the form 



The 



oo . 



(13) 



where Bi are distance independent amplitudes which cannot be obtained using our method and the correlation lengths 
are given by 



1 



1 

27 



log 



1 + e -m/T ^ 

j=0 



(14) 



The leading term of the expansion ( 13 1 is obtained considering r = in Eq. ( 12 ) and k£, satisfying 1 + e Vi ^ k o^ T = 0, 
closest to the real axis. To our knowledge, the asymptotic expansion ( [13] ) is new in the literature (the authors of 
[531 ES] did not consider the c ase of the field- field correlation functions) . Extensive numerical studies and the low- 
tcmpcrature analysis (see Sec. VIA) show that J? (l/^ rf) [«,•]) > and J? (l/^* ) [ f, »]) > for all «»(&) and Vi(k). In 
Sec. VIA we will also show thaTlTU) and (13) agree with the TLL/CFT predictions and other exact results. 



III. THE XXZ SPIN CHAIN 

The asymptotic expansions presented in the previous section were derived by taking a specific continuum limit 
in the equivalent expansions of the low-temperature transversal and longitudinal correlation functions of the XXZ 
spin chain. In order to obtain the asymptotic behavior of the correlators in the lattice model we will investigate the 
spectrum of the QTM. Therefore, it is useful to review the algebraic Bethe ansatz of the XXZ spin chain and the 
associated QTM. 

The integrable spin- 1/2 XXZ chain in external longitudinal magnetic field h is characterized by the following 
Hamiltonian 



H(J, A, h) = (J, A) — hS z 



(15) 



where 



i/(°)(j,A)=jx:[^)^ +i 

3=1 



.0-OOo-O'H 
u y u y 



1) 



(3) 



(16) 



We assume periodic boundary conditions and the number of lattice sites L to be even. The Ham iltonian (|15|) acts on 

i2l 



the Hilbert space % 



operators which act nontrivially only on the j-th lattice site <Jx,y,x 
matrices 



ir®W-l) 
a 2 



1) 



J > fixes the energy scale and A is the anisotropy. In Eq. ( |16[ ), <Jx}y,z are local spin 

with o x y z the Pauli 



















i 


°x = 


J1 




v o 
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and I2 the 2-by-2 unit matrix. S z commutes with H^(J,A) and, therefore, does not affect the integrability of the 
model. Also, due to the similarity transformation H{ J, A, h) — > VH{J, A, — /i)V -1 with V = Ylf=i a x \ it is sufficient 
to consider only the case of positive magnetic field. Another consequence is that the thermodynamics of the model 
does not depend on the sign of h. In this paper we are going to consider the massless regime of the XXZ spin 
chain |A| < 1, parametrized by A = cos 77 with < r\ < it and the magnetic field h smaller than the critical value 
h c = 8Jcos 2 (?//2). 



A. Algebraic Bethe Ansatz solution 



The Hamiltonian ( 15 ) was solved by Yang and Yang pj3"Hi?5] with the help of the coordinate Bethe ansatz. From 
the point of view of ABA [55] , which provides an alternative method of solving (15), the XXZ spin-chain is the 
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fundamental spin model associated with the trigonometric R-matrix defined by 



R(A,/x) = 



where 



/Rn(A,M) 
Ru(A,M) 
Rii(A,M) 

VRii(A,^) 



b(X,fi) 



Rx2(A,/i) 
Ri2(A,M) 
R?2(A,M) 
R?a(A,/i) 



R2i(A,M) 

R|?(a,m) 

Rli(A,M) 

R2l(A,M) 



introduction of the L-operators acting on 

2 



sinh(A — /i) 
sinh(A — /i + iff) 

it 

n = c 2 



R^(A )M )\ 

R|(a,m) 

R22(A,M) 
R22(A,M)7 



c(A,/x) 




(17) 



sinh(i?7) 



sinh(A — fj, + if?) 



(18) 



Let us show how we can obtain and solve the XXZ Hamiltonian ( 15 ) in the ABA formalism. The first step is the 

^2\<g>£ 



L,(A,0) = Yl KlK^^e^l, L,(A,0)eEnd((C 2 )^ 



(L+l) 



(19) 



a,b,ai ,bi — 1 



where is the canonical basis in 



i.e., e ab ' = e ab <g> If ^ and e%> = I 2 ® If ^ <8> e ah ® If^"^ with 
e a b the 2-by-2 matrices with all the elements zero except the one at the intersection of the a-th row and 6-th column 
which is equal to one. The additional C 2 space on which the L-operators act in addition to the Hilbert space of the 
spin chain is called the auxiliary space and, for practical purposes, it is useful to present them as 2-by-2 matrices with 
entries which are operators acting on H. Using (19) and (|17| we obtain 



Li(A,0) 



ej£ + b(X, 0)e 2 § c(A,0)e 2 { } 
c(A,0)4> 6(A,0)e«+e« 



(20) 



where is now the canonical basis in Ti. The monodromy matrix defined as 



T(A) = L i (A,0)U_ 1 (A,0)...L 1 (A,0), T(A) = ( ^ ^j) 



provides a representation of the Yang-Baxter algebra 

R(A, M )[T(A) ® T( M )] = |T(a*) ® T(A)]R(A, /x) . 



(21) 



(22) 



with R bi b2 (A,^) = (PR) bi 1 b2 2 (A,/i) = R^, 1 (A, and P the permutation matrix Pa®b = b®a for a, b £ C 2 . In Eq. (22 1, 



which signalizes the integrability of the model, T(A) (g> T(/i) should be understood as the usual tensor product between 
two square matrices of dimension 2x2 with operator valued entries as can be seen on the r.h.s. of ( pTj ). Finally, the 
transfer matrix is defined as the trace of the monodromy matrix in the auxiliary space 



t(A)=tr T(A) = A(A)+.D(A). 
In Appendix [A] it is shown that the XXZ Hamiltonian (151 can be obtained as 

H(J,A,h) = 2Jsinh(^)t- 1 (0)t'(0) - hS z . 



(23) 



(24) 



The eigenvectors and eigenvalues of the transfer matrix t(A) = A(X) + D(X) can be obtained with the help of ABA 
if we can find a pseudovacuum which is an eigenvector of A(X) and D(X) and the action of T(A) on it is triangular. 
Using the explicit expression of the L-operators in the auxiliary space (|20|) we can see that 



satisfies the requirements of a pseudovacuum and 



L times 



T(A)|0) = 



f A(X)\n) B(X)\Q) 
I C(A)|Q) D(X)\fl) 



|fi) B(x)\n) 
o (b(\,o)) L \n) 
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The operator B(X) can be interpreted as a creation operator and we will look for the eigenvectors of t(A) of the form 
|{Aj}™ =1 ) = B(Xi) ■ ■ ■ B(\ n )\H,}. In Appendix [b| we show how the eigenvalues of t(A) can be derived with the result 



-(Ai{A i} ^)=n siiih(A " Aj -" iJ7) 



sinh(A — Xj 



sinh(A 



provided that the {Aj}™ =1 parameters satisfy the Bethe equations 



i m l0\ 1 11 



sinh(A + irj/2) J 



sinh(A — Xj + irj) 
sinh(A — Xj) 



(25) 



sinh(Aj — irj/2) 
sinh(Aj + irj/2) 



n 



sinh(Aj — A s — irj) 
sinh(Aj — X s + irj) 



(26) 



In order to make contact with the results in the literature in Eqs. (25 1 and (26) we have performed the transformation 
A — » A — irj/2 and Xj — » Xj — irj/2 , j = 1, • • • , n. This also means that ( 24 ) should now be understood as H(J, A, h) = 
2 J sinh(i?7)t- 1 (it) /2)x! {irj /2) - hS z . The only thing that remains in order to obtain the spectrum of the Hamiltonian ilbl 
is to quantify the action of the magnetic operator. S z commutes with the transfer matrix therefore they share the same 



eigenvectors. Using S z = 1/2 J2j=i CT * = V 2 Xw=i( 



■ e~22 ) with e ~j the canonical basis in 



2 )® L and the explicit 
The action of the magnetic operator on the 



expression for the B(X) operator from p0| we find [S Z ,B(X)] = —B(X) 
eigenvectors of the transfer matrix is then given by S z \{Xj}™ =1 ) = S z B(Xi) ■ ■ ■ B(X n )\fl) = (L/2—n)B(Xi) ■ ■ ■ B(X n )\Q) 
where we have used S z \tt) = L/2\Vl). Finally, the energy spectrum of the XXZ spin chain in magnetic field is 



E({X}) = J2^)-h 



L 



eo(A) = 



2Jsinh 2 (z7y) 



sinh(A + irj/2) sinh(A - irj/2) 



(27) 



B. Ground-state properties 



The ground state of the XXZ spin chain at finite magnetization is constructed essentially in the same way as in the 
case of the Lieb-Liniger model. This means that the (half) integers rrij which appear in the logarithmic form of the 



Bethe equations ( 26 ) 



Lpo(Xj) - e ( x 3 - A fe) = 27rm j . 3 = !> • • • . " (28) 
fc=i 



fill all the possible values in the symmetric interval — (n— l)/2 < rrij < (n — l)/2. In Eq. (28) we have introduced the 
bare momentum po{X) and the scattering phase 0(X) 

Po(X) = iln(^^4), 0{\)=iTn( ^ + *l ), (29) 
wv ' \smh(ir)/2- X) J v ; \smh(ir] - X) J v ; 

where the branches of the logarithm are specified by the conditions linr\_ s . 00 po(X) = tt — ij and limA->oo ^(A) = n — 2r] . 
The ground state is characterized by real Bethe roots Xj which are contained in the interval [—q, q] called the Fermi 
zone. If we call every down spin a particle, then the thermodynamic limit is characterized by L — > oo , n — > oo with 
constant density of particles D = lim^ .n^oo n/L. In the thermodynamic limit the Bethe roots fill densely the interval 
[— q, q] and we can introduce the spectral density of particles p(A) which satisfies the following integral equation 

PW + f f K{X - M )p( M ) dj, = ^(A) , K(X) = 6>(X) = . Sm(2 V [ , . (30) 

2tt J_ q 2n smh(A + irj) smh(A — irj) 

The average density of particles is then D = p(X)dX from which the Fermi boundary q can be obtained. 

In the presence of a magnetic field the magnetization of the ground state is no longer fixed, it depends on the 
magnitude of h. In this case the boundary of the Fermi zone q can be defined by the requirement that the energy of a 
hole at the Fermi boundary should be zero, eo(±q) — 0, where the dressed energy £o(A) satisfies the integral equation 

£o(A) + J" K(X - n)eo(ji) dfi = h- 2Jp' (X) sinr/ = e (A) . (31) 

It can be shown that in the massless phase (|A| < 1) considered in this paper and h smaller than the critical magnetic 



field h c , Eq. (31 1 has a unique solution. When the magnetic field is vanishing the Fermi boundary goes to infinity. 



An important role in the analysis performed in Sec. [v] is played by the dressed charge Z{X) defined by the following 
integral equation 

ZW + ~ f K(X - n)Z{p) dn = 1 , Z(±q) = Z, (32) 
and the resolvent of the operator I + which satisfies 

R(\ ft) + [ K{\ - v)R(y, n) dv = ±-K{\ - M ) . (33) 

ZTT J_ q Z7T 

We will also make use of the dressed phase F(A|/i) defined by 

F(\\ti) + ^1" if (A - v)F(v\n) dv = ~0(A - M ) , (34) 
and is connected with the dressed charge via 

Z(X) = l + F(X\q)-F(X\-q), i = 1 + F(q\q) + F(q\ - q) . (35) 



A proof of the identities ( 35 ) can be found in |5SJ HH] , 



C. Continuum limit of the XXZ spin chain 

The fact that the Hamiltonian of the one-dimensional Bose gas can be obtained performing a certain continuum limit 
in the Hamiltonian of the XXZ spin chain was discovered long time ago [M] . In [53] it was shown that the Yang- Yang 
thermodynamics, Q, ((81), of the ID Bose gas can be obtained by performing the same limit in the thermodynamics 
of the lattice model derived using the QTM formalism. This is to be expected if we take into account that both 
models are integrable and that the BAEs and the energy spectrum of the Bose gas can be obtained from the BAEs 
and energy spectrum of the XXZ spin-chain in the continuum limit. Moreover, the authors of [63], see also 98 , 
derived multiple integral representations for the correlation functions of the Bose gas from equivalent expressions for 
the XXZ spin chain. In this paper we will employ a similar technique to derive the large-distance asymptotic behavior 
of temperature dependent Green's functions in the Bose gas from equivalent results for the XXZ spin chain. 

The XXZ spin chain is characterized by five parameters: lattice constant <5, number of lattice sites L, strength of 
the interaction J. anisotropy A = cost/ and magnetic field h. The Bose gas is characterized by four parameters: mass 
of the particles m = 1/2, physical length I, coupling strength c and chemical potential fi. First, we will show how we 
can obtain the BAEs of the Bose gas Q from (26 1. 



Let e — > be a small parameter. The desired continuum limit is obtained considering rj = tt — e, S — > like C(e 2 ), 
L even, L — > oo like C(l/e 2 ) with LS = I and c = e 2 /S. Performing this limit together with the reparametrization of 



the Bethe roots Xj = Skj/e in (26) we find 



cosh(-fcj +i§)\ j\ sinh(-fcj 



n 



cosh( f kj-q) J fi. sinh( f kj - f k s - ie) 

I 

1 + iffy V _ TT sinh (f % - f k s + it) 



l-i%ki) I sinh(^fc., - ^k s 



aft 
aft 



which are exactly the BAEs for the Bose gas Performing the same limits in eo(A), see (27 1, we find 



e (A) -> 2J6 2 k 2 - j^2Je 2 + ^e 4 - hj + 0{e 6 ) . (36) 

In order to obtain the energy spectrum Q of the Bose gas from (27) (we neglect the zero point energy hL/2), we need 
to consider J -> oo like C(l/e 4 ) , h -> oo like C(l/e 2 ) with 2J5 2 = 1 and // = (2 Je 2 + j e 4 - h) finite. This means that 
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by performing the thermodynamic limit followed by the continuum limit in the canonical partition function of the 
XXZ spin chain (modulo the zero point energy) we obtain the grandcanonical partition function of the Lieb-Liniger 
model 



Z 



XXZ 



(M) 



lim > i 

L->oo ^ 



-f3E({\}) 



Znls{^,/3) = lim VV 

/ — VfYl ' ■ 



HE({k}) 



(37) 



{k} 



In the following sections we will use a slightly modified scaling limit compared with the one presented before and 
utilized in [BBl. Eq. (37) can also be obtained if we consider J — 1/2, the continuum model at inverse temperature 
/3 related to the inverse temperature of the lattice model via j3 — f3/S 2 , and h — > like C(e 2 ) such that fi = 
(e 2 /6 2 + e 4 /(4<5 2 ) - h/S 2 ) is finite. Then 



/3e (A) -> 



2JS z k 



2;„2 



2Je 2 + ^e 4 - h 



P(k 2 - M ) = pe (k) . 



This shows that the thermodynamic properties and the correlation functions of the Bose gas at any temperature can 
be obtained from the thermodynamic properties and correlation functions of the XXZ spin chain at low-temperature 
and vanishing magnetic field. In the next sections we will use this continuum limit, summarized in Table [TJ to derive 
the correlation lengths of the Bose gas from the low-temperature spectrum of the XXZ-QTM. 



TABLE I. Parameters for the XXZ spin chain and the one-dimensional Bose gas. 



XXZ spin chain 


One-dimensional Bose gas 


lattice constant 8 — 0(e 2 ) 
number of lattice sites L = 0(1/ e 2 ) 
interaction strength J = 1/2 
magnetic field h — 0(e 2 ) 
inverse temperature /3 
anisotropy A = cos rj — e 2 /2 — 1 


particle mass m = 1/2 
physical length / = L8 
repulsion strength c = e 2 /8 

chemical potential fi = [e 2 /8 2 + e i /(4S 2 ) -h/S 2 ) 
inverse temperature /3 = (38 2 



IV. THE LOW-TEMPERATURE SPECTRUM OF THE XXZ SPIN CHAIN QUANTUM TRANSFER 

MATRIX 

In this section we are going to investigate the low-temperature spectrum of the XXZ spin chain QTM [SSI G2] • 
A short review of the relevant facts about the QTM can be found in Appendix [C] The QTM is important for two 
reasons: first, the largest eigenvalue, which we will denote by Ao(A), completely characterizes the thermodynamics of 
the system via 

/(fe,T) = -ilogAo(0), 

with f(h,T) the free energy per lattice site. In general, the largest eigenvalue of the QTM can be expressed in 
terms of some finite number of auxiliary functions satisfying non- linear integral equations |99j . This is a very efficient 
thermodynamic description for the model contrasting with the Thermodynamic Bethe Ansatz (TBA) 100 approach 
which relies on the string hypothesis and provides an infinite number of NLIEs. The second reason is given by the 
fact that the correlation lengths of various Green's functions can be obtained as ratios of the largest and next-largest 
eigenvalues of the QTM 1101] , Appendix [Cj This is a consequence of the finite gap between the largest eigenvalue 
and the rest of the spectrum of the QTM. In the next sections we are going to study the low-temperature spectrum 
of the QTM in order to obtain the asymptotic expansion of the longitudinal and transversal correlation functions in 
the XXZ spin chain. Performing the continuum limit presented in Sec. |III C| we are going to arrive at the results 
presented in Sec. II A 
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We introduce two types of L-operators L,(A, -u') , Lj(u', A) e End ((C 2 )« , ( JV + 1 ') denned as 

a,b,ai : b\ = l a,b,ai,fei— 1 



where vf = iu , u — — 2Jsinry-^, TV is the Trotter number and e^, the canonical basis in (£. 2 )®( N+1 '> . The monodromy 
matrix of the QTM, which provides another representation of the Yang-Baxter algebra (22), is given by 

t qtm (a) = Ljv(Aj _ U ')L 7V _ 1 (?/. A) ■ • • L 2 (A, -u'^u', A) . 
Using the explicit expression of the L-operators in the auxiliary space 



_ / e« + b(X, -u>)e% c(A, - U ')41 } ^ i r ,/ u _ f + &(«', A)e« c(u', A)e« 



it is easy to see that 



to- J • ; J • J , (39) 



N factors 

satisfies the conditions of a pseudovacuum for the monodromy matrix of the QTM and 

jQ™m\n) ( AQTM W\n) s QTM (A)|o) \ _ / (b(u',x)) N/2 in) s«™(A)|fi) \ U()] 

WW y C Q™(X)\Q) DQ™(X)\Q) J [ (6(A, -u')) N/2 p) J ' [ ' 



The presence of the magnetic field in the Hamiltonian (|15j) is easily taken into account by the following transformation 
of the monodromy matrix 

T«™(A)->T«™(A)( ^ mV (41) 
V e 2 J 

The quantum transfer matrix t ( ^™(A) is defined as the t rac e in the auxiliary space of the monodromy matrix 
t c ?™(A) = tr T QTM (A). The existence of the pseudovacuum (39 1 and the fact that T QTM (A) provides a representation 
of the Yang-Baxter algebra ensures that the eigenvalues of the QTM can be obtained using the ABA. Making use of 
(|40| , (Ell) and the results of Appendix IbJ the solutions of the eigenvalue equation 



t°™(A)|{A}) = tQ™(X)BQ™(X 1 ) ■ ■ ■ BQ™(X p )\n) = A(A)|{A}) , 

are given by 

A(A) = b(u', Xfl\^ f[ S ^^M + b(X, - M r /2 e" W f[ Sin - n hfA A - 3 A + r ) ' (42) 
provided that the parameters {Aj}j =1 satisfy the Bethe equations 

b{Xj,-u') J ax S mh(Aj - X k - irj) 

As it is shown in Appendix [C] (see also [551 HOlj ) the asymptotic expansion of the longitudinal correlation function 
is given by 

(a( 1) cr( m+1) ) T = const + £ , m -> oo , (44) 
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where A; are unknown amplitudes, l/^ d) = log(A (0)/A^' t; (0)) and the sum is over all the next-largest eigenvalues 
Aj Pft, ^(0) in the N/2 sector. We remind the reader that an eigenvalue of the QTM is said to be in the M sector if 



iph). 



p = M in Eqs. (42) and (43). The asymptotic expansion of the transversal correlation function is given by 



/ (1) (m+l)\ \ D 

{a\'al_ ') T = } J Bj e 



(45) 



where Bi are unknown amplitudes, l/^ (s) = log(A (0)/A,f } (0)) , a { £ = {a { J ] ±ia [ y j) )/2 and the sum is over all the 

next-largest eigenvalues A^ (0) in the N/2 — 1 sector. 

The QTM method can also be utilized to investigate the generating functional for the a z correlators, i.e., 
( e l»'En=i e 22 }) Tj from which the longitudinal correlation function can be obtained via (o* o* m )y = (2D^ n d 2 — 



AD m d v + l)(e^^™=i e 22 })t| v =o with D m the lattice derivative defined as D m a m = 
(in)neN- ln this case the asymptotic expansion is 



for any sequence 



} )t = 



00 , 



(46) 



In (46) (see Appendixjc]) the sum is over all the eigenvalues of t$™(0) = A Q ™ (0) + e v D Q ™{0) in the N/2 sector 
denoted by A,[ v) (0) and l/^ v) = log(A (0)/A| v) (0)). 



A. Non-linear integral equations for the largest eigenvalue 



Deriving NLIEs for the largest and next-largest eigenvalues of the QTM requires a different method from the one 
utilized in the computation of the ground state and low-lying excitations of the transfer matrix in the massless regime. 
This is due to the fact that in the Trotter limit, N — > oo, the distribution of Bethe roots (43) in the complex plane 



presents an accumulation point at the origin and isolated solutions which makes it impossible to introduce the "density 
of roots" like in the case of the ground state of the transfer matrix. Fortunately, the Bethe roots appear only in some 
strips of the complex plane which allows for the introduction of some auxiliary functions which satisfy functional 
equations. Thanks to fundamental properties of the gross distribution of {Aj } P =1 the auxiliary functions enjoy certain 
analyticity properties which allow to transform the functional equations in terms of non-linear integral equations. 
Eventually the eigenvalues of the QTM can be expressed in terms of these auxiliary functions. This method, which 
was developed in [53 [58l 11021 1103] . will be our main tool in investigating the spectrum of the QTM. 
The largest eigenvalue of the QTM lies in the N/2 sector. We will employ the following notations 



4>±W 



■ v, f\ _i_ ■ \ \ N / 2 N / 2 

Smh(A±m) ) , ,(A) = n^(A-A J ), 
sin 77 / ±L 

i / j= i 



which allows to express the eigenvalues of the QTM ( 42 ) as 

A (A) 



</>-(A) q(X~irj) eh 



0+(A) q(X + irj) e _ m 



0_(A-i»/) <?(A) <MA-Mt/) g(A) 

Below, the NLIE and integral expression for the largest eigenvalue will be determined following 



(47) 



Integral equation for the auxiliary function 



An extremely important role in the following will be played by the auxiliary function o(A), which is periodic of 
period in and defined by 



<MA) (j)+{X + iri) q(X - irj) 



We note that the Bethe equations (43) can be rewritten as o(A 
a(A) 



T has 3A^/2 solutions in a period strip, of which, only N/2 are given by the Bethe roots {Aj}.-^. The 



(48) 



,N/2. However, the equation 

N/2 
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irj/2 



ITT 



in/2 



ir]/2 



277/2 



c 



o o 



o 



o 



o 



o 



o 



x--y 



j 

iu 



o 



o 



o 



o O 
o O 
o 



iri- 



FIG. 1. Typical distribution of Bethe roots (•) and holes (o) in the strip \^sX\ < r\ , r\ £ [0, 7r/2) characterizing the largest 
eigenvalue of the QTM at low-temperatures. All the other roots and holes can be obtained using the in periodicity. The 
contour C surrounds all the Bethe roots and the pole of the auxiliary function a(A) at iu. 

additional N solutions are called holes and we will denote them by {Xj }jL\- A typical distribution of Bethe roots 
and holes for -q € (0, 7t/2) and low temperatures is presented in Fig. [I] Let C be a rectangular contour with positive 
orientation, centered at the origin, extending to infinity, with the upper (lower) edges parallel to the real axis through 
±i(rj — e)/2, with e — > 0. It is important to note that this contour is independent of the Trotter number and the 
following considerations are valid for all TV. Inside the contour, the function 1 + a(A) has N/2 zeros at the Bethe roots 
and a pole of order N/2 at iu, which means that log(l + a(A)) has no winding number around the contour allowing 
us to define (A is located outside of C) 



/(A) 



1 

2wi 



dX 



(log sinh(A - /«)) log(l + a(fi))dfi 



1 

2n~i 



log sinh(A — fx) 



q'(M) 

1 + a(n) 



dfi . 



(49) 
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For the evaluation of the r.h.s. of Eq. (49) we will use the following theorem: 

Theorem 1. \101$ Let C be a contour in the complex plane, and let g(X) be a function analytic and non-zero inside 
and on C. Let (f>(X) be another function which is analytic inside and on C except at a finite number of points; let the 
zeros of <fi(X) in the interior of C be a\, a,2, ■ ■ ■ and let their degrees of multiplicity be r±, r^, • ■ ■ ; and let its poles in the 
interior of C be 61,62)' ' ' an d let their degrees of multiplicity be sx, s%, • • • . Then 



iezeros tepoles 



obtaining 



N/2 

/( A ) = X! lo S sinn ( A _ A J") _ y logsinh(A - iu) = logg(A) - log^_(A) - — log sin?/. (50) 



Eq. ( |50| p rovides an integral representation for q(X) in terms of log(l + a(A)). Taking the logarithm in Eq. (48) and 
using ( |50| we can derive 

which is a nonlinear integral equation of convolution type for the auxiliary function a(A), valid for all N. Using 

lim log [ J = -2J/3sinh(i77)cothA, (51) 



JV^oo °\<f>_(\) 

the Trotter limit, N — > 00, can be performed with the final result 

„ 2Jsinh 2 (i?7) 1 f sin(2r7) 

log a(A) = -ph - ■ ' r - / -^ttt ■ 7 U \ ^ log 1 + a ^)W ■ (52) 

smn(A + in) sum A In J c smn(A — fx + vr\) smn(A — \x — ir\) 



Eq. (52) was obtained under the assumption that r\ <= (0, 7r/2). It is also valid for 77 G (n/2,n) but in this case C is 
a rectangular contour centered at zero, extending to infinity, with the upper (lower) edges parallel to the real axis 
through ±i(n — rj — e)/2 with e —> 0. 



Integral expression for the largest eigenvalue 

It remains to obtain an integral expression for the largest eigenvalue Ao(A) in terms of the auxiliary function ct (A) . 
Consider n <E (0, tt/2) for which the distribution of roots and holes is presented in Fig. [I] First, we note that Eq. (47) 
can be rewritten as 



A Q (A) 



P(A) 



4>-{\- in)<f> + (X + in)q(X) ' 



(53) 



with p(X) — 0_(A)^_|_(A + irj)q(X — in)e (ih l 2 + </> + (A)0_(A — in)q(X + ir\)e ( ih l 2 . The function p(X) is quasi-periodic 
p(X + in) = (-l) 3N / 2 p{X) and lim A ^ 00 p(A)/(sinh A) 3JV / 2 = (e^ 1 / 2 + e-^ h ' 2 )/{s\nr}) N . The zeros of the function p(X) 
in a strip of width in are the solutions of the equations a(A) = — 1 yielding 



N/2 



N 



pw = c n sinh ( A - a j) n sinh ( A - a ? } ) ' 



(54) 



where {Aj}^f 2 and {X^}^ =1 are the Bethe roots and holes, respectively, and c is a constant. Defining q^(X) 
\~\f =1 sinh(A — Xj ) and using (54) to replace p(A) in Eq. (53) we obtain 



A (A) 



Jh) 



(A) 



0_(A - ir))<f>+(X + in) ' 



(55) 
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which provides an alternative expression for the largest eigenvalue in terms of the holes and not Bethe roots. Below, 
we will show how an integral representation of log q( •'(A) in terms of the auxiliary function a(A) can be calculated. 

We consider a new rectangular contour with positive orientation C (see Fig. [T]) extending to infinity, with the upper 
(lower) edges parallel to the real axis through i(rj — e)/2 and —i(rj — e)/2 + iir with e — > 0. The lower edge of the 
contour C at i{rj — e)/2 coincides with the upper edge of C but has opposite orientation. Now we can prove the 
following identity 



d(X - At) . a , . dfi = , d(X-ji) = ^- logsinh(A - /i) . 



c+c 



1 + o(m) 



dX 



(56) 



First, we notice that the contributions of the two contours parallel to the real axis through i(jj — e)/2 cancel each 
other due to the opposite orientation. Then it can be easily verified using the definition of a(A) (48 1 that the functions 
appearing in (56 1 are all periodic of period in, which means that the upper and lower edges of C+C do not contribute 
to the integral. Finally, the contributions of the sides parallel to the imaginary axis are also zero as it can be seen 
from 



lim d(X — /i) 

Sft/i— >+oo 



_^L = m 1 lim 1 Um 

l + a(fi) o(/x) 1 + K/*->±°° 1 + a _1 (M) 1 + e^' »/i->±«> a(|u) 



0. 



Consider A close to the real axis. Then making use of (56 1 we find 



1 

2ttI 



c 



d(A-M) aV \ ^ = - — 



d{X-^)-^-d^. 
1 + o(m) 



(57) 



The r.h.s. of (57) can be calculated using Theorem [T] by taking into account that the function 1 + a(A) has inside the 

N/2 simple poles at {Aj}^f^ + irj and a pole of order 



contour C, N zeros at the holes {A^}^ or {A^}- 
N/2 at — iu — irj + in with the result 



N 



1 

2?ri 



1 + a(/x) 



TV 



AT/2 



^d(A-Af')-E d ( A 



AT 

\j - irj) - —d(X 



iu - 



irj) 



(58) 



Using again Theorem [I] and the fact that the function 1 + a(X) has inside the contour C, N/2 zeros at the Bethe roots 



{Xj}f=i and a pole of order N/2 at iu we find 



1 

2~n~i 



d(X — ji 



ir i) TTT7^ rf ^ = 51 d ( A ~ A i ~ ir l) 



N 



1 + a(/i) 



d(A 



(59) 



The integral representation for logg^^A) is obtained by taking the difference of Eqs. (58) and (59), integrating by 
parts, and then integrating w.r.t. A with the result 



1 

2ni 



[d(X - fi) -d(\-(j,- irj)] log(l + a(n))dfi = - log q {h) (A) + log (c/5 + (A + irj)<j>- (A - ir?)) + c 



(60) 



where c is a constant of integration. Making use of this integral representation in Eq. ( 55 ) the largest eigenvalue of 
the QTM can be written as 



logA (A) =c- 



1 f sinh(i?7) 
2ni Jc sinh(A — fi — irj) sinh(A — fi) 



log(l + a{ji))dji. 



The constant of integration is computed using the behavior of the involved functions at infinity. Performing the change 
of variables z = X — ji in the integral, and using lim^oo log A (A) = log(e^/ 2 + e~^ h ^ 2 ) and lim^oo log(l + a(A)) = 
log(l + e~@ h ), we find that the constant of integration is (3h/2. The final result for the largest eigenvalue of the QTM 
evaluated at is 



logA (0) = 



1 

2n 



sm-q 

sinh(^i + ir/) sinh(^) 



log(l + a(ji))dji. 



(61) 



The integral expression (61), which was obtained for rj £ (0,7r/2), is also valid for rj £ (n/2,n) but, as in the case of 
the NLIE for the auxiliary function ( 52 ) , the contour C should be replaced by a similar rectangular contour with the 
upper (lower) edges parallel to the real axis through ±i(n — rj — e)/2 with e — > 0. 
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3. Final form of the integral equations 



The NLIE (52) and integral expression (61) are in fact correct for all temperatures [99]. This is due to the fact that, 
even at high temperatures, the Bcthc roots are contained in the strip |SA| < (77 — e)/2 for 77 g (0,7r/2) (or the strip 
|3A| < (ir — r] — e) /2 for r\ G (tt/2, 7r)), which means that the reasoning of the previous sections is still valid, producing 
the same results. However, in order to obtain the thermodynamic properties and correlation lengths of the Bose gas 
we will be interested only in the low-temperature limit in which some simplifications of (52) and (61 1 appear. 

Let us consider 77 G (tt/2,tt). In this case, the upper edge of the contour C, which we will call C + , is a parallel line 
to the real axis through i(ir — rj — e)/2 (for the following discussion the presence of the e term is irrelevant). Then for 
A G C+ , A = x + i(TT — r\)/2 with x real, the driving term on the r.h.s. of (52) is negative and equal to 



~Ph- 



cosh(a; 



2 J sin 77 
irj/2) cosh(a; 



i V /2) 



which means that in the low-temperature limit /3 — > 00 , (/i, J > 0), the contribution of the upper part of the contour 
is negligible and we can restrict the free argument A and the integration variable to the lower part of the contour. We 
can shift this line to the line parallel to the real axis through —177/2 without crossing any poles of the driving term 
obtaining 



log a(\-ir)/2) = ~Ph~P- 



2Jsinh 2 (i77) 



1 



sin(2??) 



sinh(A + 777/2) sinh(A — 777/2) 2n J K sinh(A — ji + 777) sinh(A — fj, — irf) 



log(l+o(/^-7?7 /2))dfi . 

Applying a similar reasoning to the integral expression of the largest eigenvalue (61 ), after the shift at —777/2, we find 

Ph 1 f sin?/ 



logA (0) 



2 2tt 7 r sinh(/_i + 777/2) sinh(/_7 — 777/2) 



log(l + a(/j - ir//2))d(i . 



(63) 



Let us introduce the function e(A) satisfying e e ( A '/ T = a(A — 777/2) where T = 1/P is the temperature. Then noticing 
that the driving term in Eq. ([62]) is the magnon energy (271, and using (29) and (30), the NLIE for the auxiliary 



(64a) 
(64b) 



function and the integral expression for the largest eigenvalue at low-temperatures can be written as 

e(A) = e (A) + A j K (A - fx) log (l + e~^' T ) dia , 



logAo(O) = 



h 1 

2T + 2^ 



P ' (X) log (l + e - £ ^/ T ) dX. 



These equations are very similar to the Yang- Yang equation for the excitation energy and the grandcanonical potential 
of the Bose gas PQ . Following [S3] , in Sec. VI we will show that the Yang- Yang thermodynamics can be obtained from 
Eqs. (64) if we perform the scaling limit presented in Sec IIIC Even though Eqs. (64) were obtained for 77 6 (7r/2,7r), 
they are also valid for 77 G (0,7r/2), which can be proved by using appropriate contour manipulations. As these are 
beyond the scope of this paper, we confine ourselves to assuming the validity. This assumption will be verified in 
Sec.|V|where we will show that they reproduce the TLL/CFT predictions for the free energy and asymptotic behavior 
of the correlation functions. 



B. Integral equations for the next-largest eigenvalues in the N/2 sector 

Computing the correlation lengths for the Green's function (<t* ai m+1 ^)T requires the derivation of integral equa- 
tions for the next-largest eigenvalues of the QTM in the N/2 sector. This means that, as in the case of the largest 
eigenvalue, we will have N/2 Bethe roots and N holes. In the previous section we have derived the NLIE for the 
auxiliary function and the integral expression for the largest eigenvalue making use of the fact that the Bethe roots 
were located in a relevant strip (modulo the periodicity) of the complex plane which was independent of the Trotter 
number and temperature. In the case of the next-largest eigenvalues in the N/2 sector at low-temperatures, some of 
the Bethe roots are found outside of this strip and an equal number of holes are inside the strips. We will employ 
the same method used in Sec. |IV A| but modified in such a way that these Bethe roots and holes are properly taken 
into account. The calculations are presented in Appendix [D] At low-temperatures, the next-largest eigenvalues of the 
QTM in the N/2 sector are given by 

logAf^O) = A +*5>(Ai) - <5>(AJ) + ^ jU(A)log (l + e -MA)/x) dX ; (65) 
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with the auxiliary functions Uj(A) satisfying the NLIEs 

Ui(X) = eo(A) - iTY,9(\ - \+) + iT^ 0{X - Xj) + £ jf if (A - /*) log (l + e 



-udp,)/T 



dfi . 



(66) 



The parameters {A^~}J =1 and {A ? }J = i, which belong to the upper, resp., lower half-plane are fixed by the constraint 
1 + e~ Ui( - X i )/T = 0. In Eqs. (65) and (661, r can take the values 1, 2 • • • . The subscript i enumerates the sets of 
parameters {A^}J =1 satisfying the constraint 1 + e - "*^ ^ T = 0. The eigenvalues with the next-largest magnitudes 
are obtained for r — 1 and the A parameters closest to the real axis. 



C. Integral equations for the next-largest eigenvalues in the N/2 — 1 sector 



The next-largest eigenvalues in the N/2 — 1 sector are relevant for the computation of the correlation lengths of 
the Green's function )t- The eigenvalues in this sector are characterized by p = N/2 — 1 in Eqs. (|42|) and 



(43 1. Some of the features encountered in the study of the largest and N/2 sector eigenvalues are also present m this 
case. At low-temperatures, the next- largest eigenvalue in this sector has N/2 — 1 Bethe roots and a hole in a certain 
strip of the complex plane. Eigenvalues with decreasing magnitude are obtained by moving pairs of Bethe roots/holes 
outside/inside the strip. This means that it is sufficient to obtain integral equations for the case with one hole in the 
strip, the equations for the other eigenvalues are obtained by adding extra driving terms of the type encountered in 
Eqs. ( |65| and (66). The necessary calculations are presented in Appendix [E| At low-temperatures, the next-largest 
eigenvalues in the N/2 — 1 sector have the integral representation 



logAHO) 



2T 



— iir + i 



ij2p (A+)-ij2p (Xj) + ±- f ^(A)log(l + e-^ A )/ T ) dX. 

3=0 3=1 77 JR 



(67) 



with the auxiliary functions Mi (A) satisfying the NLIEs 

Vi(X) = e (X)±iTrT-iTy r 6>(A-A+)+iTV 9(X-XJ) + — f K (A - p) log f 1 + e ~ v ^ )/T ) dp . (68) 

~k ~, J 2tt J r V ) 



3=0 



3=1 



The 2r + 1 parameters {A^}^ =0 and {A Yj=i which belong to the upper, resp., lower half-plane are fixed by the 



when Aq~ is in the first (second) quadrant of the complex plane 5RAJ > , 3A|j~ > (3?Aj < , 3AJ > 0). Compared 



constraint 1 + e~ Vi ^ x i ^ T = 0. On the r.h.s. of Eq. (68) the plus (minus) sign in front of the iirT term is considered 
e fi 

with Eqs. (|65 1 and 
xiditions 

of Eq. (|68|. The next-largest eigenvalue in this sector is obtained for r = and A J closest to the real 



(|66j) , the integral expressions and NLIEs for the next-largest eigenvalues in the N/2 — 1 sector 
present an addition al, u npaired parameter Aj, an — in factor in the r.h.s of Eq. (67) and an additional factor of UnT 
in the r.h.s. 



axis. 



V. COMPARISON WITH THE TLL/CFT PREDICTIONS 



In Sec. IV A 3 we have derived an integral expression (64), for the largest eigenvalue of the QTM in terms of an 
auxiliary function which obeys a NLIE very similar to the Yang- Yang equation (fsl) . Eq. ( 64 ) and the similar integral 



representations for the next-largest eigenvalues (65) and (67), are valid only at low-temperatures and, in the course 



of the derivation, we have made some assumptions which were not fully justified for some values of the anisotropy. 
Here, we will show that our results are in perfect agreement with the predictions of the Tomonaga-Luttinger liquid 
[1051 1106] and Conformal Field Theory [4"51 11071fll2j . confirming the validity of our assumptions. 



A. Low-temperature behavior of the free energy 



At low-temperatures CFT predicts |112j that the free energy per lattice site scales like 

7rT 2 r h f q 

f(h,T) = e (h)~—^ + O(T 3 ) 7 €o (h) = --+ e (X)p(X)dX, 
6v F 2 J_ a 



(69) 
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where eo(h) is the density of the ground state energy, c is the conformal charge, (not to be confused with the coupling 
constant of the Bose gas), which is equal to one in the case of the XXZ spin chain, and vf is the Fermi velocity defined 
as 



vf 



e'o(q) 



(70) 



with e' (q) the derivative of the dressed energy (31 ) evaluated at the Ferm i boundary q. 

Let us show that the free energy per lattice site obtained from Eq. (64b), via f(h,T) = — TlogAo(O), satisfies ( |69| . 
Performing an analysis similar to the o ne in Appendix A of [T] or Chap. I of |38| it can be shown that for r\ € (0, it) 
and h < h c = 8Jcos 2 (ry/2), the NLIE (64a I for the auxiliary function e(A) has two zeros on the real axis which we 



will denote by ±q(T). Also e(X) is negative on (—q(T),q(T)) and positive outside of this interval. Let us denote 
lim-T^o l(T) = q and limj^o £ W = ^oW- Then using 



lim Tlog fl 



-e(A)/T 




A € (-q,q) , 

A e (—oo, —q) U (q, +oo) , 



(71) 



we find that in the low-temperature limit the NLIE (64a) transforms in the linear equation for the dressed energy 
(31 1. The equation for the dressed energy has a unique solution for 77 € (0, tt) which means that lhriT->o e (X) = £o(A) 



and limx^o<z(r) = q. In order to show that Eq. ( 64b[ ) gives the correct free energy satisfying (69) we need a more 
accurate estimation of integrals containing the factor log (l + e _e ^ A ^ T ). In the following we are going to assume that 
for low temperatures the auxiliary function has the expansion [7S1 11131 1114j 



e(A) = e (A) + e x (A)T + e 2 (A)T 2 + 0(T 3 ) . 



(72) 



Then it can be shown (see Appendix A of [55] or |115j ) that, for any function /(A), bounded on the real axis and 
differcntiable in the vicinity of ±q we have : 



lim T 

T->0 



/(A) log (l + e-^' T ) d\=- T f(X)e(X)d\ 



TV <m + /(-,)) + T ^m + T2£2 t q K { - q) + o(T3) 



6e' (q) 



24(9) 



(73) 



We should mention that a more compact, but maybe not as transparent, method of investigating the low-temperature 
limit of the QTM spectrum was employed in |101j . All the results derived here and in the following sections can also 
be obtained utilizing the results of the aforementioned paper. Using ( 73 ) and substituting ( 72 ) in the equation for 
the auxiliary function (64a) we obtain 



2 2 
E^ £fe (A) + i-][> fe f K ( x -V) 

k=0 k=0 "? 



ek{p) dfj, = e (A) 



^ (K(X -q) + K(X + q)) -<!)- T ^ 7 ^K{X + q) + 0(T 3 ) 



l2e' {qY Ane'^q) 
Equating terms of the same order in temperature we find 



iwe' (q) 



e 1 (A)=0, e 2 (X) 



6e' (q) 



(R(X,q)+R(X,-q)), 



(74) 



with the resolvent i?(A,/x) defined in (33). The low-temperature expansion of the free energy per lattice site is 
calculated using the asymptotic formula ( 73 ) in Eq. ( 64b ) with the result 



f(h,T) = - 



h 



1 

2n 



3 ttT 2 
p' (X)(e (X) + £2 (A)T 2 ) dX - ^r^P'M + 0(T 3 ) 



(75) 



where we have used the fact that p' (X) is even, p' (— q) = p'o(q)- Using Eq. (74) and the identity 



p' (X)R(X,±q) dX=p' (q)-2wp(q), 



(76) 
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(for a proof, see Appendix [f]), (75) takes the form 

h 1 



f(h,T) = - 



2 2tt 



(77) 



Using the identity ( F3 ) and the definition of the Fermi velocity ( 70 1 , it is easy to see that this expression is identical 
with (1691). 



B. Low-temperature asymptotic behavior of the longitudinal correlation 

The Tomonaga-Luttinger liquid theory and CFT predict the following asymptotic behavior of the longitudinal 
correlation function at low-temperatures (we consider only the leading-order of the oscillatory terms) [48] 



{TZ/v F f 



2 sinh (irTm /i>f) 



J2 A 



i e 



ttT/vf 



s'mh(TTTm/vF) 



2l z Z 2 



m 



oo , 



(78) 



with the Fermi momentum defined by kp = ftD and Ai are coefficients that do not depend on T. The analysis of the 
QTM spectrum in Appendix [O shows that the asymptotic behavior of the correlation function can be expressed as 



r (lWm+l)\ _ 



) T = const + lie « <d) i 



oo , 



(79) 



where the sum is over all the correlation lengths, \/^ d \ui] = log(Ao(0)/A^(0)), determined as the ratio of the largest 



and next-largest eigenvalues in the N/2 sector. Using ( |64b I and (65) we obtain the following explicit expression for 
the correlation lengths 



i 

~2tt 



p' (X) log 



i(A)/T 



1 



-e(\)/T 



dA-i£>,(A+) 

3=1 



3=1 



(80) 



where the functions e(A) and lij(A) satisfy Eqs. (64a) and (66). In the rest of this section we will show that (79) is 
equivalent to (78) in the conformal limit. 

The analysis of the correlation lengths (80) in the limit T — > is very similar to the one performed by Kozlowski, 
Maillet and Slavnov for the correlation lengths of the Bose gas, and, for this reason, we are going to use some of the 
notations and terminology employed in |55j . In the following we are going to discard the subscript i for the auxiliary 
function Uj(A). We are considering an arbitrary auxiliary function u(A) satisfying Eq. (66), with 2r parameters 



{Aj~}J =1 , located in the upper (lower) half-plane, which also satisfy the constraint l + e~"^ A j = 0. The 

first observation that we are going to make is that limr-^o u(X) — £o(A). In analogy with the case of the auxiliary 
function e(A) we expect that all the solutions of the equation 1 + e~ u ( A )/ T = will collapse at ±q in the limit T — » 0. 
We will say that the solutions that collapse at q (— q) belong to the right (left) series. If we assume that u(X) has the 
expansion 



u(A) = e (A) + ui(A)T + u 2 (A)T 2 + 0(T 3 ) 



(81) 



then a formula similar to ( 73 1 can be derived as in [55 



lim T / /(A) log f 1 + e - u( - x ^ T ) dX = - I f(X)u(X)d\ 
Jr v > J-q 



+ ££m + /(-*» + T ^m + T2u2 t q K^ q) + or) . 



64(g) 



(82) 



We are going to consider that the roots {A 1 * 1 } are distributed in the following manner: r+ roots A + and roots A - 
belonging to the right series (collapse at q); r~ roots A + and roots A~ belonging to the left series (collapsing at 
—q), where and satisfy the constraints 



r v - r t = r h - r P = L 
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with I integer, satisfying — r < I < r. More explicitly, at sufficiently low temperatures we have 



{At}r =1 = {g + iTa+}£ x U {-g + iTaj;}^ , R(a*) > , 



(83a) 
(83b) 



{AJ}5 =1 = {g - iT^}^ U {-g - /Vlllr , > , 

where a£ and ^ satisfy 1 + e -«(±9+iTa±)/T = i + e -«(±g-iT/?±)/T = q The leading Tay i or coefficients can be 
parameterized by a set of integers, p k and via 

~), u(±q-iTp±) = T 2mT{s±-^). 



u(±q + iTa±) = ±2niT( P ± 



Using the expansion (81), e(±q) = and e' (—q) = —s' (q) wc find 



± _ 2tt / ± 1\ ,ui(±g) 



4(g) 



2?r 



4(g) 



M±g) 

T 4(g) 



(84) 



Substituting the parametrization ( 83 1 in Eq. ( 66 ) and expanding the driving terms up to the second order in T we 
find 



(A) = eo(A) + ^ / K(\ - fi) log (l + e -«W/ r ) dfx + 9l (X)T + g 2 (X)T 2 + 0(T 3 ) 



(85) 



with 



and 



<)2 



gi (\) = -il(6(\-q)-6(\ + q)) 



(A) = —K(X q) £ <4 + E - *(A + g) E ^ + E & 



.fc=i fe=i 



, fe=i 



fe=i 



We can now use (82 1 in (851 obtaining a system of equations for the unknown functions ui(A) and U2(A). For u\{X) 
it reads 



1 

2tt 



if (A - fi) ui(fi) dfi = -il(9(X -q)- B{\ + q)) . 



Comparison with the integral equation for the dressed phase (34) shows that iti(A) = — 2nil(F(X\q) — F(X\ — q)). 

>f1 

«?(«) 



Using the first identity in (35) we can obtain an expression in terms of the dressed charge iti(A) = 2iril(l — Z(X)). 
The equation for U2(A) is 



«2(A) 



2tt 



K(X - n) u 2 (fx) dfi = 52(A) + K(X - q) 



124(g) 47r4(g) 



+ K{X + q) 



* "l(-g) 
124(g) 47r4(g) 



with the solution 



u 2 (A) = -i?(A,g) 



1 fix 



k=l k=l 



-i?(A,-g) 



1 /vr 



fc=i 



fc=l 



24(g) V 3 

Using (73) and (82) in (80) and expanding the j>o(A ) terms up to the first order in T we obtain 
1 ^u\(q)p(q) 



24(g) V 3 



•«?(-?) 



m F - + 27TTp(q) ( e ^ + E ^ + E ^ + + E ^ I + °( y2 



(86) 



. fe=i 



fe=l fe=i 



fe=i 



In deriving (86) we have used also the identity f_ p' Q (X)Z(X) dX = J p(X) dX and (FI ). Finally, using (84) we find 



-2ilki 



2-kT 

Vf 



(87) 



fc=i k=i fe=i fe=i 



The second term in the expansion (78) is obtained for r = 1,1 = and pf = Pl =1 (or s x = s x = 1). The next 
leading terms are obtained for r = I , I = 1, 2, • • ■ and the integers and (or p^" and st) taking values from 1 to I. 
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C. Low-temperature asymptotic behavior of the transversal correlation 



In the case of the transversal correlation TLL and CFT predict the following asymptotic behavior at low- 
temperatures [48] 



l e 



ttT/vf 



smk(-KTm/vF) 



m — > oo . 



with Bi coefficients that do not depend on T. The analysis of the correlation functions in the framework of the QTM 
in Appendix [O showed that 



{<?+ o_ )t 



m — > oo . 



(89) 



where the sum is over all the correlation lengths = log(Ao(0)/Ai^(0)), determined as the ratio of the largest 

and next-largest eigenvalues in the N/2 — 1 sector. Using (64b I and (67) we obtain the following explicit expression 
for the correlation lengths (we neglect the in term which produces an (^l) m factor) 



i 

'2tt 



Po(A)log 



1 + e -ViW/T^ 
I + e -e(X)/T 



' 3=0 3=1 



(90) 



where the functions e(A) and Vi(X) satisfy Eqs. (64a) and (68). It is sufficient to consider the conformal limit of the 
following correlation length 



1 



1 

'2tt 



Po(A)log 



l + e -*(A)/T\ 

i + e -eW/T ) dX-i P0 (X+). 



(91) 



with v(X) satisfying the equation (68) with r = 0. The behavior of the correlation length (90) will be obtained by 
summing the contributions of (19 1 1) and ( 80 1 derived in the previous section. 



First, we notice that limr->.o v(X) — £o(X), therefore, we are going to consider the following expansion 

v(X) = e (A) + wi(A)T + v 2 (X)T 2 + 0(T 3 ) , (92) 



with Vi(X) and V2(X) unknown functions. Also (82 ) is still valid if we replace u(X) with w(A). We are going to consider 
that An" is located in the first quadrant (this means that we are going to have a plus sign in front of the iirT term in 
Eq. (68)). The same result is obtained if we consider Aj in the second quadrant. Then at sufficiently low-temperatures 
we have 



q + ia^T, a£ 



2tt 



Po 



+ i 



Mi) 



(93) 



with po an integer parameterizing the leading Taylor coefficient ■ Using this parametrization in Eq. ( 68 1 , and 
expanding to the second order in T we find 



(A) = e (A) + ^ / K(X — ii) log (l + e -^)/ T ) dfi + [m - i9(X - q)]T - K(X - q)a+T 2 + 0(T 3 ) 



(94) 



The integral equations satisfied by fi(A) and v 2 (X) are obtained by substituting (82), modified for the v(X) function, 
in (94). For vj(X) it reads 



»i(A) 



2tt 



K(X — fj,) Ui(/z) e?/i = iir — i9(X — q) . 



Comparison with the integral equation for the dressed phase ( 34 ) and dressed charge ( 32 ) shows that v\ (A) 
inZ(X) — 2inF(X\q)). Using the first identity in (35) the solution to this equation can be rewritten as i>i(A) 



47r(l - F(X\q) - F(X\ - q)). The equation for v 2 (X) is 



v 2 (X) 



1 

2^ 



K(X - (i) u 2 (m) dpi = K(X - q) -aj + 



12e' (q) 4ne' Q (q) 



+ K(X + q) 



* «j(-g) * 

12e' (q) 4ne' (q), 
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with the solution 

v 2 (X) = -R(X,q) 



-2ira2 



1 



1 



•«?(-?) 



Using (73 1 and (82) in (91) and expanding the pq{Xq) term up to the first order in T we obtain 

1 



0(2") 



In d erivi ng ( |95| ) we have also used the fact that F{X\q) + F(X\ — q) is an odd function of A, (F(— X\ — n) 
and (Fl ). The final result follows from (93) and the use of the second identity in (35) 



1 



£(«)[«] 



2nT 



+ P0-1 +0{T 



The asymptotic behavior of the correlation length (|90|) is then 

2-kT 



1 



1 



(95) 
(96) 



- , = —lilkw , , 
£( S >M v F \ AZ 2 



(97) 



fc=0 



fc=i 



fc=i 



fc=i 



The leading term of the expansion (88) is obtained for r = I = and po = lj which means that Aq is the closest 
solution to the real axis of 1 + e~ Vi ( x o ^ T — located in the first or second quadrant. The next leading terms in the 
same expansion are obtained for r = Z , / = 1, 2, ■ ■ ■ ,po = 1, and the integers and sjT (or p^ and s~t) taking values 
from 1 to /. 



VI. CONTINUUM LIMIT 



In the previous sections we have obtained NLIEs and integral representations for the auxiliary functions and 
eigenvalues of the QTM in the N/2 and N/2 — 1 sector, valid at low-temperatures. In 63J it was shown that by 
performing the continuum limit presented in Sec. |III C| the Yang- Yang thermodynamics of the one-dimensional Bose 
gas can be obtained from the largest eigenvalue of the QTM. The next natural step is to perform the same limit in 
the equations for the next largest eigenvalues obtaining the spectrum of what we can call the "continuum" quantum 
transfer matrix. The ratio of the largest and next-largest eigenvalues of this "continuum" QTM will give the correlation 
lengths of the density-density and field-field correlation functions of the Bose gas. The correspondence between the 
correlation functions in the two models is presented in Table [TTJ It should be noted that the results obtained for the 
Bose gas are valid at all temperatures and are not restricted to low-temperatures as in the case of similar results 
obtained for the XXZ spin chain. 



TABLE II. Correspondence in the continuum limit between the correlation functions of the XXZ spin chain and the one- 
dimensional Bose gas. 



XXZ spin chain 


One-dimensional Bose gas 


(4 1 Vi m+1) ) T 


O'0»0i(o)>T 

(* t (x)*(0)) T 

/ e ffo JO') 



We start by showing how we can obtain Q, ([8]) from (64). Performing the continuum limit presented in Sec. IIIC 
(this also includes the reparametrization A = 5kje ) we obtain 



p (X)^Sk, p' (X)^e, 9(X) 



-6{k) 



9'(A) = K(X) 



2c 



5 k 2 + c 2 



\K{k). 



and eo(X)/T — > eo(k)/T, where T is the temperature in the Bose gas. Using these results and denoting by e(k)/T the 
continuum limit of e(X)/T we see that Eq. (64a) transforms into the Yang- Yang equation (18]). The grandcanonical 
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potential of the Bose gas per unit length </>(//, T) is related to the free energy of the Heisenberg spin chain per lattice 
constant through the relati on (j)((i,T) = (f(h,T) + h/2)/S 3 . Using f(h : T) = -TlogA (0) it is easy to see that (R 
can be obtained from (64b) in the continuum limit. 



Wc define the eigenvalues of the "continuum" QTM by 

logA(0) = i(logA(0)- A 

where on the r.h.s. of this relation the continuum limit of logA(O) is understood. For the eigenvalues in the A^/2 
sector we obtain 



logAf ft) (0) 



3 = 1 



with Ui(k) satisfying (|9| and for the eigenvalues in the A^/2 
is irrelevant in the continuum limit) 



1 

2^ 



log 1 



-Ui{k)/T 



dk . 



(98) 



1 sector we find (the iw term on the r.h.s. of Eq. (67) 



i°gA< a) (o) = i£*+-X>7 + i / log ( 



l + e 



i(k)/T 



dk , 



(99) 



with Vi{k) satisfying ( 12 ). The correlation lengths of the density-density correlation function, (j(x)j(0))t, are obtained 
as ratios of the largest "continuum" eigenvalue, log A (0) = ^- / R log (l -I- e~ e ( k ^ T ) dk, and the next-largest "contin- 
uum" eigenvalues in the N/2 sector justifying (11). In the case of the field-field correlation function, (x)$(0))t, 
the correlation lengths are obtained using the next-largest eigenvalues in the N/2 — 1 sector with the result ( [14] ). The 
case of the generating functional is treated in Appendix [G] 



A. Checking the results 



The asymptotic expansions (10) and (13) which are valid for all temperatures should reproduce the TLL/CFT 
results [49 1 150 ] : 



(TZ/v F ) 



2$mh ji (TrTx/v F ) 



Ai e 2ixlkf 



I el 



ttT/vf 



smh{-KT x / v f) 



2l 2 Z 



X — > 00 . 



(100) 



<*t (a; )* ( o)} T = £s 



irT/vp 
smh(TTTx/vp) 



(101) 



in the T — > limit. In Eqs. (100) and (101), v f is th e Fermi velocity defined in (H4), kp = n D is the Fermi momentum 



uated at q (see H2 ). The agreement with the conformal results is proved in Appendix 



and Z is the dressed charge eva 

El 

In the impenetrable limit, c — > oo, the leading term of the asymptotic expansion for the field-field correlation function 
was computed by Its, Izergin and Korepin by solving an associated Riemann-Hilbert problem |32| . Chap. XVI. of 
[38], This gives us another opportu nity to check the validity of our resu lts by comparison with another exact result. 
The leading term in the expansion ( Jl3[ ) is obtained when r = in Eq. (12). We will consider that fcg~ is in the first 
quadrant, SJfcg > , > 0. Taking into account that 

lim K(k) = lim 9{k) = 



the equations (12) for the auxiliary function v{k) (we drop the subscript i) and dressed energy (|8j) become 

v(k) = k 2 — [i + inT , £(fe) = k 2 — jU . 



(102) 



The asymptotic behavior depends on the sign of the chemical potential. We will consider first the case of negative 
chemical potential. In this case the solutio n of t he equation 1 + e~ v ^ k( > ^ T which is closest to the real axis and located 



in the first quadrant is fcg = i\J Using (102) and this value for , the correlation length (14) can be rewritten as 

1 1 f , / e {k2 



1 

27 



log 



1 



o{k*-n)/T _ 1 



dk 



/i < 
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which is precisely the result obtained in [33] for negative chemical potential. In the case of positive chemical potential 
we have k^ = ^fji and the correlation length is 



1 



i M 



coinciding with the result derived in |32j . 



1 

27 
1 

2tt 



log 
log 



/ e (fc 2 - M )/T + l 
e (fc2- M )/T _ I 

e (fc 2 - M )/T , l 



dk 



3 (fe 2 - M )/T _ X 



e?fc . 



VII. CONCLUSIONS 



Using the spectrum of the XXZ spin chain QTM and a specific continuum limit we have derived the asymptotic 
expansions of the temperature dependent density-density and field-field correlation functions in the interacting one- 
dimensional Bose gas. As a by-product we have also obtained similar expansions, valid at low-temperatures, for the 
longitudinal and transversal correlation functions in the XXZ spin chain. One could naturally expect that similar 
results can be derived in the case of the spinorial ID Bose gas |117j which can be obtained as the continuum limit of 
the U g (sl(3)) Perk-Schultz spin chain |118U119] . This subject will be deferred to a future publication. 
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Appendix A: Derivation of the XXZ Hamiltonian from the transfer matrix 



In this appendix we are going to show th at t he logarithmic derivative of the transfer matrix is identical with 



H(°\j, A) (modulo a constant) proving Eq. (24 1. We introduce two types of operators, which we will call R- and 
P-operators, acting on (C 2 )®( L+1 - ) and defined by 

R;, fe (A,//)= £ Kt^^/aL' Pi,k= E P tllAei% 2 , (Al) 

a\ ,fai ,02 ,&2 — 1 <3l ,bl ,&2 — 1 

where P^ 2 = 5a 1 h 2 ^a,6 1 and is the canonical basis in (C 2 )®( L+1 - ) . The L-operators are a subset of the R-operators 
Lj(A,0) = Ro.j(A,0) and P^fc acts on an arbitrary vector in (C 2 )®( L+1 ' by permuting the j-th and fc-th component. 
Some useful properties of the (permutation) P-operators are 

Pj > iP*j = Pfc,iPi,J = p wP*,i. Pfc,i = p i,fc. p L=I, = e%Pj, k . (A2) 

First, we are going to calculate t(0) = tr o T(0). Using the fact that R(0,0) = P, Lf(0, 0) = Pqj and applying 
successively the first identity of ( |A2[ ) we find 

T(0) =P 0j lPo,l-i"-Po,2Po,h 

=Pl,l-iPl,l-2 • • ■ Pl,2Pl,iPo,l i 



P| .Pj..!---P/. I /P./.. 

Taking into account that tr Pq,l = I we obtain t(0) = Pi i 2P2,3 • ■ • Pl-i,l with t _1 (0) = Pl-x,l ' 1 1 P2.3P1.2- Fo r the 



computation of the derivative of the transfer matrix t'(0) = troT'(O), we will also need the last identity of (A2| and 
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the fact that Pj,kd\Ri,m(\ 0) = <9 A Rz im (A, 0)Pj,fc if j, k j^l,m. We find 

L 

T'(0) = Y, P o,i ' ' ' P<M+i 9 A Ro,i(A, 0)| A=0 P 0j< _i • • • Po.i , 



i=l 
L 



^ Po : L • • • Po,i+lPo,i-l 9xRi-l,i(A, 0)| A=0 • • • Po,l , 
i=l 
Z, 

Po : l • • • Po,i+iPo,i-i • • • Po,i <9 A R,_i i i(A, 0)| A=0 , 

L 

Po,lPl,L ' ' ' Pl,i+lPl,i-l ' ' ' Pi, 2 9 A Ri_i ! i(A, 0)| A=0 , 

i=l 
L 

Po,lPl,2 ■ • ■ Pi-l,i+lPi+l,i+2 ' ' ' Pl~1,L <3\Ri-l,j(A, 0)| A=0 , 



i=l 

therefore, t'(0) = tr o Tg(0) = Yh=\ p i,2 • • • Pi-i,i+iP*+i,i+2 ' ' ' p l-i,l d\Ri-i,i{\, 0)| A=0 . Collecting everything we 
obtain 

L L 

t(0)-V(0) =^Pi- M flARi-i,i(A,0)| A=0 =^9AR i _i, i (A ! 0)| A=0 , 

t=l i=l 

where periodic boundary conditions Ro,i = Rl,i ar e understood. Using the explicit expression for the R-matrix it is 
easy to see that 2Jsinh(i?7) 9 A R i _i ) j(A, 0)| A _ Q is equal to J[ax &x +&y <Ty + A(o-i I_1 Vi^ — 1)] proving Eq. (24 1. 



Appendix B: Algebraic Bethe ansatz for the generalized XXZ spin chain 



Let us consider a general monodromy matrix which is intertwined by the R-matrix (17 1 

R(A, m)[T(A) ® T( M )] = [TO*) <8 T(A)]R(A, M ) , (Bl) 
and assume the existence of a pseudovacuum such that 

jrfvv|o\ = ( A (m) B(X)\n) \( a(A)|Q) B(A)|fi) \ 

1 ;| ' \c{X)\n) D(A)jn) / \ d(A)|n> y ' 

The monodromy matrices of the XXZ spin chain and quantum transfer matrix are particular cases of this model with 
o(A) = 1, d(A) = (fo(A,0)) L and a(A) = &(it', A) A '/ 2 e^' i / 2 , d(A) = 6(A, -u') Ar/2 e~' 9 ' 1 / 2 respectively. We are interested 
in finding the eigenvalues of the transfer matrix t(A) = A(X) + D(X). Interpreting the -B(A) operator as a creation 
operator we are going to look for eigenvectors of the type |{Aj}^ =1 ) = B(X 1 ) ■ ■ ■ B(X p )\Q,) satisfying the eigenvalue 
equation 

t(A)|{A i }5 =1 )=r(A|{A J -}? =1 )|{A i }? =1 ). 

This eigenvalue equation will impose a set of equations on the {Aj}^ =1 parameters which we will call Bethe equations. 
Before we derive these equations, we need the commutation relations between the operator valued entries of the 



monodromy matrix. Using the explicit matrix representation of Eq. (Bl ) 

fx o\ 

c(A,/x) 6(A,m) 
b(X,n) c(X,fi) 
\ 1 / 



( A(X)A(n) A(X)B(n) B{X)A{ix) B(\)B(p)\ 

A{X)C{p) A(X)D(p) B(X)C(n) B(X)D( f i) 

C{X)A(p) C(A)B( M ) D{X)A{h) D{X)B{h) 

V C(A)C( M ) C(X)D(p) D{X)C{fx) D{X)D{p) j 

( A{n)A{X) A{fi)B(X) B(n)A(X) B{^)B(X) \ ( 1 \ 



A(p)C(X) A(p)D(X) B(ji)C(X) B(p)D(X) 
C((i)A(X) C( M )S(A) D(p)A{X) D(n)B(X) 
\C(ji)C(X) C(fi)D(X) D( f i)C(X) D{p)D{X) ) 



c(A,/i) b{X,fi) 
b(X,fi) c(X,fi) 

V 1 / 
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the following commutation relations can be obtained from the matrix elements (1,4), (1,3) and (2,4) 

B(X)B(n) = B(f,)B(\) , 

A(X)B(n) = f( P , X)B(fi)A(X) - h{fjL, X)B(X)A(n) , 
D(X)B(n) = /(A, n)B(fj,)D(X) - h(X, ^)B{X)D{p) , 



where we have introduced 



/(a,aO 



b(X,n) ' 



h(X,fi) = 



c(X,n) 
b(X,fx) 



(B2) 
(B3) 
(B4) 



(B5) 



Now we can evaluate the action of the transfer matrix t(A) = A(X) + D(X) on the prospective eigenvector 
B(Xi) ■ ■ ■ B(Xp)\tl). Acting with A(X) on B(X\) ■ ■ ■ B(X P )\Q) and using the commutation relation (B3) we obtain 
2 P terms of which only one has the desired form 



Y[f{X j ,X)a(X)B(X 1 )---B(X p )\Q). 



This term is obtained by using p times only the first term on the r.h.s. of (B3 1. The other 2 P — 1 terms which should 
cancel in order to satisfy the eigenvalue equation can be written as 



Mj(X, {A}) J B(A 1 ) • • • B(Xj) ■ ■ ■ B(X P )B(X)\Q) , 

3=1 

where the hat means that the corresponding operator is missing. At first impression obtaining the coefficients 
Mj(X, {A}) may seem a daunting task, but, fortunately, the commutativity of the B(X) operators simplifies the 
matter c onsi derably. The Mi (A, {A}) coefficient is easy to compute: it requires to use first the second term of the 
r.h.s. of (B3) and then the first term p — 1 times with the result 



-a(Ai)/i(Ai,A)n/(A fe ,Ai; 



k=2 



Now using the commutativity of the B(X) operators we argue that the Mj(X, {A}) coefficient can be obtained from 
Mi (A, {A}) substituting Ai with Xj obtaining 



M,(A,{A}) = -a(Xj)h(Xj, A) JJ f(Xk,Xj) ■ 
Therefore, the action of ^4(A) on B(X) ■ ■ ■ B(X p )\ft) can be written as 



(B6) 



A(X)B(X 1 ) ■ ■ ■ B(X p )\n) = a(X) ]J f(X j ,X)B(X 1 ) ■ ■ ■ B(X p )p) + J^M^X, {X})B(Xx) ■ ■ ■ B(Xj) ■ ■ ■ B(X p )B(X)p) , 

(B7) 

with Mj(X, {A}) given by (B6). In an analogous fashion the action of D(X) on -B(Ai) • • • B(X p )\Cl) is computed as 

p p 
D(X)B(X 1 ) ■ ■ ■ B(X P )\Q) = d(X) n /(A, A,)B(Ai) ■ • • B(X P )\Q) + £ N,(X, {A})5(Ai) • ■ • B(Xj) ■ ■ ■ B(X P )B(X)\Q) , 

3=1 3=1 

(B8) 

with 



A^AJA}) = -h(X,Xj) n f(Xj,Xk)d(Xj 



(B9) 



Then Eqs. (B7) and (B8) show that B(Ai) • • • B(X m )\ty is an eigenvector of t(A) = .4(A) + D(A) with eigenvalue 



r(A|{A}P =1 ) = o(A) I] /(Ai, A) + d(A) f[ /(A, A j 



2G 



if the parameters {\j}^ =1 satisfy the equations Mj(X, {A}) + Nj(X, {A}) = , j = 1, • • ■ ,p or, equivalently 

= -rr /(Aj.Afc) 
d(A,) /(Afe, Aj) ' ' ' 

The following remark is in order: Strictly speaking our treatment is not completely rigorous since we have implicitly 
assumed the linear independence of vectors of the form -B(Ai) • • • B(X p )\fl). The interested reader can find the rigorous 
solution of this problem in Chap. XII of |120j . 

Appendix C: The XXZ Spin Chain Quantum Transfer Matrix 

In this paper the asymptotic behavior of the temperature dependent correlation functions is investigated using the 
quantum transfer matrix method. In this appendix we review the principal features of the method in order to make 
the paper self-contained. The considerations below follow closely the presentation in |60j . 

We introduce N auxiliary spaces denoted by 1, • • • ,N and two types of monodromy matrices using R-operators (see 
Appendix [A| : 

Typel T J (A) = R Ji£ (A,/i)-.-R J)1 (A,M), (Gla) 
Type 2 T (A: R ; ,://. A: • • • R, ,://. A: . (Clb) 

Using the results of Appendix |A"| we have 

t(0) = trjTj(O) = Pi, 2 P 2 , 3 • • • Pl-i,l , 

and 

L 

t'(0) = trjU(0) = 5^Pi J a-"P<-i,<+iPi+i,<+2---Pi-i ) L 0aR*-i,<(A,Q)| a=o , 

i=l 

which implies that 

L L 

1(0)^(0) =$>i-i,i flARi-i,i(A,0)| A=0 =^9 A R i -i )i (A ! 0)| A=0 , 

i=l i=l 

2 J smh(ir]) ' 

under periodic boundary conditions Rq.i = Rl,i- The last relation shows that 

In a similar fashion we can show that 

t(0) = trjTj(O) = P L , L _! ■ ■ • P 3 . 2 P 2 ,i = t _1 (0) , 

and 

L 

t'(0) = trjTj(O) = d x Ri,i-i(0, A)| A=0 Plx-i • ■ ■ P<+2,i+iP<+M-i ' ' ■ p 2,i , 
i=i 

implying that 

L L 

t , (0)t" 1 (0) = Y, »xRi,i- 1 (0,A)| A=o P M _i = -J2 ^Ri,i-i(A,0)| A=o , 



1=1 



g""(J,A) 
2Jsinh(i?7) 
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where we have used Rk,jPj,k = Rj,kRj,k = Rjk and the fact that 9^Rj ) j_i(A, 0)U =0 = — d\Ri^\{Q, A)U_ . Again, 
from the last relation we see that the following expansion is valid 

^ = cxp (- A 2TSi + °( A2 )) ?(0) - (C3) 

Employing a special form of the Trotter-Suzuki formula with v! = — 2 Jsinh(i7y)-^ 

/ 1 \ N 

lim (t- 1 (0)t(u / )) = I™ ( 1 + -^(—(3H(°\j, A) + 0(1/N)) ) = e -0" (o> W A > , (C4) 



we find 

AT/2 



lim p N x 



N— >oo 



= (Ef-tOtK))^ = (^l + ^(-/3^ )(J I A) + O(l/7V))j =e -^°V,A)_ (C5) 



This formula is the starting point for the introduction of the QTM. Using Eqs. (CI) we have 

PN,L =^J t ...,N=l { J n(- U ') T N^i( U ') ' ■ ■ J 2(- U ') J r( u ')} » 
= tr T,.- ,N-1 {T]v(-«') T 7v=t( m ') • ■ ■ T2(-w') T T<y )} , 
=tr T,.- ,AT=i {^atO^ -^O^lvCMi -«') • • • R i,-i,jv(M' -"ORi.lvCM) -«') 

R|n>' m) R^t >2 ( u » • • • r£=i («', m) R^t , £ («', m) 



Rl,2(M) -W') R 2,2(M= ~ u ') ■ ■ ■ R Z,-l,2(Mi " u ')Rl,2(M, -"') 

R X j l( u/ '/ i ) R t :2 ( n ''/ x ) ' ' ' R T : L-i( u ')M)Rx i (u',/i)} 
H { R !,1v(m, -wOR^t^K^) ■ ■ • Rl,2(M, -u')^^,!!) 



M=0 



"* r l,--- ,N- 



R L,n(^ ~ U ') R F3T,l( U '' ^) " • • R L,2(M, - u ') R T iL ( U '' ^)} 

^{T?™(0)--Tg™(0)}, (C6) 

where (R* 1 )^ = R^ and we introduced the monodromy matrix of the QTM 

Tf™(\) = R^(X,~u')R t ^ lj (u\\)---R jr2 (X,-u')Rl j (u',X). (C7) 



The QTM is defined as t Q ™(\) = tr j Tf™(\). The usefulness of Eq. (fC6|) will become obvious if we notice that 



Zxxz(P)=ti 1 ,.., L e-^'^ = lim tn,... , L p N , L 

N— >-oo 

= lim tr 1 ,..., L tr T ... ] ,(Tr / (0)---Tr / (0)) 

iV— !-oo ' ' \ / 

oo 

= lim tr T ... jj (tQ™(0)) L = E A « (°) > ( C8 ) 



where the last sum is over all the eigenvalues of the QTM. In the case of the XXZ spin chain QTM the following 
assumptions are true: 

• The limits L — > oo and N — > oo are interchangeable [S51 1121] . 

• The largest eigenvalue of the QTM which is denoted by A (A) is real, positive, nondegenerate and separated by 
a gap from the next- largest eigenvalues in the Trotter limit N — > oo. 
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Therefore, using ( C8 ) the free energy per lattice site of the XXZ spin chain is 

y logZ xxz (P) logAo(O) f . 

fW = * JSoo JL = J- ■ (C9) 

This result shows that the thermodynamic behavior of the lattice model is completely determined by the largest 
eigenvalue of the QTM evaluated at providing an elegant solution to the almost impossible problem of summing 
over all the eigenvalues of the Hamiltonian. 

Until now we have considered the case without magnetic field. The presence of the magnetic field operator S z = 
\ Y^tf=i a< z^ m the Hamiltonian |l5| can be easily taken into account in the QTM formalism. The first observation is 
that 

/ e Ph/2 o \ 

[R(A !M ),e®e] = o, e=(* e -eh/2)> ( C1 °) 

with the obvious consequence that O <g> is a spectral parameter free solution of the Yang-Baxter algebra with 
R-matrix (fiB, i.e., R(A, (i)(Q ® 6) = (6 ® 0)R(A, fi). Therefore 



lim p N L e pl%bz = e 

JV->oo 



/3hS z _ p -P(H^(J,A)-hS z ) ^ 

(j?™(o)e 1 ---T'i™(o)&L) , (en) 



which shows that the presence of the magnetic field term in the Hamiltonian is taken into account by the following 
transformation 

jQTM{x) ^ jQTM{x J \ ( 

\ e 2 / 



1. Correlation functions within the QTM approach 

The QTM method provides considerable simplifications in the treatment of temperature dependent correlation 
functions. We are interested in correlation functions of local operators of the following type 

(°i "-Ofc-^T- bm Le _ pmJAih) , (C13) 

where j, ■ ■ ■ , k € {1, • • • , L} with j < k are the local spaces on which the operators — if ® O (8 lf^ L ~^ act. 
Using ( C6 ) we have 



( , w tr^tr,.., ( T r J (0) • • • Tf™(0)O« • • • 0[% +1 



l.w^oo tri... L e-^ ff ( J ' A '' 1 ) ' 

tr T ... F {(^(Ojy-^CT^COjOi) • ■ • tr(T«™(0)O fc _ J+1 )(t^(0)) i - fc } 



= lim 



z,,jv-i-oo tri...x, e~P H( - J ' A ' h ) 

Um (W |tr(TQ TAf (0)00 ■ ' ' tr(TQ™(0)O fc _ J+1 )|vl/ ) (nji 



Eq. ( C14 1 was obtained assuming that the QTM is diagonalizable and has "normalized" eigenvectors denoted by |^ n }- 
The eigenvector j^o) corresponds to the largest eigenvalue Aq(0). This assumption is valid in the case of the XXZ spin 
chain [j but in the case of other lattice models the rigorous mathematical proof is lacking. From Appendix [b] we know 



2 F. Gohmann, private communication. 
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that the eigenvectors of the QTM are of the form |*> = B Q ™ ' {X x ) ■ ■ ■ B Q ™(X p )\fl). We will say that an eigenvector 
\9) is in the p sector if in the previous expression there are p "creation" operators B^™(X) (C"^™(A) can be 
interpreted as a "destruction" operator as a result of C^™(A)|fi) = 0). We will also assume that the eigenvectors in 
different sectors are orthogonal. Now we can consider some particular cases of Eq. ( [CT4| . 



Transversal correlation function. First, we will consider the case 0\ = cr_ , Oi,-~ ,k-j = ^2 ,Ok-j+i = er + . Using 



tr(T Q ™(0)(7_) = BQ™(0) and tr(T Q ™(0)a+) = C®™{0), Eq. (C14| becomes 



U) w = (*q\BQ™(0)(AQ™(0) + g^W^g^ 
[ - + }T n^oo A^' +1 (0) ' 

Taking into account that j^o) is in the N/2 sector and the interpretation of the C^™(X) as a "destruction" operator 
we are going to assume that the vector = C^™ (0)\^>q) and its adjoint (^'| = (^fo\B^™ (0) can be expanded as 

i£^r — 1 sector itz^f — 1 sector 



Using these expansions in (C15) we find 

£ £ -tt — 1 sector 



A (0) 

(fc-3) 



— 1 sector 



J2 B t e 4 a) , fc»j, (C16) 



where i?i are unknown constant coefficients, = log(Ao(0)/A^ s (0)), with A^ (0) the eigenvalues of the QTM in 

the N/2 - 1 sector. 

Longitudinal correlation function. In this case 0\ — a z ,02— k-j = I2,0fc-j+i = o~ z and tr(T M (0)<r z ) = 
AO™(0) - DQ™(0). We obtain 

w (Jk) = (tt o |(AQ™(0) - DQ™(0))(AQ™(0) + ^ M (0)) fc -^ 1 (AO^(0) - pQ™(0))\*o) 

° z >T n^L Ag" J+1 (0) ' 1 j 

The vector |*') = (A®™ (0) - £>Q™(0))|* ) and its adjoint (*'| = (* |(^ Q ™(0) - £> Q ™(0)) can be expanded as 

|*')= ]T am, <*'|= J] c,(^|, 



which means that 



zG4f sector z£ 4f sector 



-iG 4f sector 



A (o) 



const + B i e ' > fc >J- ( C18 ) 

i£ sector ,2^0 



In the last line of Eq. (C18l the constant term is the contribution of the largest eigenvalue (which lies in the A^/2 
sector) and the sum is over all the eigenvalues in the N/2 sector denoted by A^ h) (0) with l/^f = log(A (0)/A^ ?i) (0)). 



Comparing with the conformal result (78) we can identify the constant with ((Tz^)t- 

Generating functional. The generating functional for the a z correlation functions is obtained for 0\ = 0% ■ ■ ■ = 

Ok-i+i = I ^ ^ I = e ve22 . For j = l,k = m we obtain 
0+ \0 e* ) J 

(e w^A?} )T= lim m^m+^^m (C19) 

V ' N^oo A™(0) ' V ' 
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In this case compared with ( |C15| or( |C17[ ) in ( |C19[ ) appears A®™ (0) + e ¥ ' J D < 5™(0) instead of t Q ™(0) = 
A ( 3™(0) + D®™(0). This does not pose a serious problem because t$™ = A c ?™(0) + e^D^™(0) is the so- 
called twisted QTM which can be easily solved by noting that the considerations of Appendix [B] applies also in the 
case of t$™ with o(A) = b(u', \)*/2 e Ph/2 ^ d ^ = b ^ _ u /)jv/2 e -/jtya+ V- Assuming that 



|5>n 



(¥>)! 



(C20) 



sector i£ sector 

where the sum is over the N/2 sector of the twisted QTM we find 

y c K'"»Y" 

VJ m»l, 

z£ ^ sector 



(C21) 



where we have denoted by A{ v) (0) the eigenvalues of the twisted QTM in the A/2 sector and 1 /£f> = log(A (0) /A^> (0)) 



(v) _ 



M, 



Appendix D: Derivation of the integral equations for the next-largest eigenvalues in the N/2 sector 

1. Integral equation for the auxiliary function 



For reasons of clarity we are going to consider first the simplest case in which only one Bethe root /hole is out- 
side/inside the relevant strip in the complex plane. The generalization to the case of r pairs is a natural extension 
of this particular example. A typical distribution of roots and holes for 7/ G (0,7r/2) and low-temperatures is pre- 
sented in Fig. pi where we have denoted by A - and A + the Bethe root, respectively, the hole outside (inside) the 
strip |3A| < 77/2. It should be emphasized that this "particle-hole" distribution is valid only at low-temperatures, at 
higher temperatures the next-largest eigenvalues in the A/2 sector are characterized by the so-called 1-string type 
and 2-string type solutions 122 . The eigenvalue and the auxiliary function a p h{X) corresponding to the distribution 
presented in Fig. [2] are described by the same formulas as in (47 1 and (48 1. It is useful to present q{\) in the following 
form 



9(A) = sinh(A - A") TJ sinh(A - Xj) , 



where {Aj}^^ 1 are the N/2 — 1 Bethe roots inside the strip |3A| < rj/2. The equation a p /j(A) + 1 = has 3A/2 

solutions, of which A/2 are the Bethe roots and N are holes denoted by {Xf^fj-^ and A+. 

We introduce the rectangular contour C p h (see Fig. k|) centered at the origin, extending to infinity with the edges 
parallel to the real axis through ±i(i] — e)/2 , e — > which presents an indentation of the upper edge such that A~ +if] 
is not in the interior of the contour. Inside the contour C p t the function 1 + a p h(X) has A/2 zeros at the Bethe roots 
{^■j}f=i 1 an d hole A + and a pole of order A/2 at iu. Therefore, the function log(l + a p h{X)) has no winding number 
around the contour (the presence of the indentation ensures that the function log(l + a p h(X)) does not have an extra 
pole at A~ + irj) allowing us to define (A is located outside the contour C p h) 



f P h{X)= — / — r (logsinh(A — /j,)) log(l + a p h(^))dfi 
which can be evaluated using Theorem [l] with the result 



1 

2iri 



log smh(A - n) 
c ph l + OphW 



N/2-1 N 
/ph(A) = E l°gsinh(A — Aj) + logsinh(A — A + ) — — logsinh(A — iu) . 
3=1 



dfi, (Dl) 



(D2) 



Eq. (D2) which can be rewritten as logg(A) = f p h( A ) + logsinh(A — A ) — logsinh(A — A + ) + log <fi-(X) + N/2 log sin ij 
providing an integral representation for log q( A) . Taking the logarithm of Eq. ( 48 1 and using this integral representation 
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FIG. 2. Typical distribution of Bethe roots (•) and holes (o) in the strip |3A| < 77,77 £ [0,7r/2) characterizing one of the 
next-largest eigenvalues of the QTM in the N/2 sector at low-temperature. All the other roots and holes can be obtained using 
the in periodicity. The indentation of the contour C v h excludes A~ + 777 and a hole located close to it. The lower edge of the 
contour C' ph is the same as the upper edge of C p h but with different orientation. The contour C p h surrounds all the Bethe roots 
except A~ the hole located at A + and the pole of the auxiliary function a p h(A) at iu. 



we find 
loga p/l (A) 



</>+{X) 0_(A + n7) 
<MA) 4>+{X + iii) 



-log 



/ sinh(A — A + iff) 
\sinh(A — A~ — irf) 



-log 



sinh(A — A + + irf) 
sinh(A — A+ — irf) 



+fph(^+iv)-f P h(^ 



Performing the Trotter limit, N — > 00, with the help of Eq. (|51|) we obtain the NLIE for the auxiliary function 

K{\- ii)\og{l + a ph {ii))dii. (D3) 



log a ph (X) = -pe (X + irj/2) + i0{\ - A+) - id{X - A") 



1 

2^ 
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Eq. (D3) was obtained assuming 77 £ (0,tt/2). It remains valid also for 7/ £ (7r/2,7r) if the contour C p ^ is replaced by 
a similar rectangular contour with the upper (lower) edges parallel to the real axis through H{n — 77 — e)/2 , e — > 
but, in this case, without the indentation. 



2. Integral expression for the next-largest eigenvalue in the N/2 sector 



The integral expression for the next-largest eigenvalue in the N/2 sector is obtained in a similar fashion as in the 
largest eigenvalue case. The starting point is, again, the representation ( 55 ) of the eigenvalue in terms of the holes 
where it is useful to denote q^(X) as 



N-l 



q {h) (X) = sinh(A - A+) Y[ sinh(A - A/ ) 



In order to obtain an integral expression for gW(A) we introduce a rectangular contour C' ph (see Fig. \2n extending 
to infinity with the edges parallel to the real axis through 1(77 — e)/2 and —7(77 — e)/2 + in. The edge at 7(77 — e)/2 
presents an indentation such that A_ + 777 is contained in the interior of C' ph and is identical with the upper edge of 
the contour C p h but with opposite orientation. Then the following identity 



L 



d(X - /i) 



*'ph(v) 



1 + Uphill) 



dfi = , 



(D4) 



can be proved in exactly the same way as its largest eigenvalue counterpart ( 56 ) . For A close to the real axis using 
( |D4| and Theorem [l] we obtain 



1 

2ni 



'N-l 



d i x - *j H) ) + d i x - A~) - d(\ - A" 

J =1 



N/2-1 
i=i 



irj) 



N 



irj) 



d{\ + iu + irj) 



(D5) 



In deriving Eq. (D5|, we have used the fact that, inside the contour C' h , the function 1 + a p / 1 (A) has: N zeros at 

A~ + in, {Xj or {A^}^ 1 + in , N/2 — 1 simple poles at {A^}^^ 1 + irj, a simple pole at A~ + ir\ and a pole of 

order N/2 at — iu — 777 + in. Using again Theorem [I] and the fact that inside the contour C p h the function 1 + a p h(X) 
has: N/2 zeros at the Bethe roots {X^f^i 1 and hole A + , and a pole of order N/2 at iu we find 



1 

2ni 



d(X — li — irj) 



u'phil 1 ) 
1 + Ophill) 



N/2-1 



d[i = d(X — Xj — ir/) + d(X — X + — imj) d(X — iu — irj) 



(D6) 



Taking the difference of Eqs. (D5| and (D6), integrating by parts, and then integrating w.r.t. A we obtain the following 
representation 



\ogq {h) {X) =log 



sinh(A - A+) 
sinh(A — A+ — 777) 
1 

2™ j Cp! 



log 



sinh(A- A~) 



sinh(A — A — 777) 
[d(X - fi) — d(X /_/ 777)] log(l + a p h{p)) dpt + c 



log (<p + (X + irj)(j>- (A - 777)) 



(D7) 



with c a constant of integration. Finally, the integral expression for the next-largest eigenvalue of the QTM in the 
N/2 sector is obtained by replacing (D7| in (551 with the result 



logA p/l (0) = ^-+log 



sinh A^ 



sinh(A+ + 777) 



- log 



sinh A 



sinh(A + 777) 



1 

2~ji 



pG(M + *»?/2)log(l + OphO*))dM. (D8) 



The constant of integration, j3h/2 1 was calculated using the behavior of the involved functions at infinity, like in the 
case of the largest eigenvalue. Eq. (D8| is also valid in the domain 77 £ (n/2,n) if the contour C p h is replaced by a 
rectangular contour, extending to infinity, with the edges parallel to the real axis through ±i(n — 77 — e)/2 , e — > 0. 
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3. Final form of the integral equations 



In Sec. IV A 3 we have shown that for rj € (7r/2,7r) the contribution of the upper contour in the NLIE (52 1 can be 
neglected as a consequence of the fact that the driving term is very large and negative. A similar reasoning can be 
applied in the case of the NLIE (D3). Restricting the free parameter A and the variable of integration to the lower 



part of the contour, and then performing a shift to the line with imaginary part —ir\/2 we find 

1 



log a ph (X - irj/2) = -j3e (X) + iO(X - A+) - iO(X - A") - 



2tt 



K(X - n) log(l + Ophiv - ir)/2))dn , 



(D9) 



where we took into account that A + — s- A + —irj/2 , A — >• A — ir\/2. The expression (D8 1 for the next-largest eigenvalue 
becomes 



logA p/l (0) = ~ + ip (X + ) - ip (X ) + ~ / p' O) log(l + a ph (n - w?/2))d/i. 



(D10) 



Introducing the function u(X) satisfying e "l A )/ T = a p h(X — ir]/2), the NLIE for the auxiliary function and the integral 
expression for the next-largest eigenvalues in the N/2 sector at low-temperatures can be written as 



u(X) = e (A) - iT9(X - A+) + iT6(X -X~) + ^- [ K(X - fj,) log f 1 + e - u ^ )/T \ dfi 

2n ,/» V / 



(Dll) 



logAp^(O) 



2T 



+ i Po (X+) - ipo(X-) + i- J^p' (ri log (l + e- u ^/ T ) d/i , 



(D12) 



where the parameters X^ 1 satisfy the constraint 1 + e _u ^ A± ^ T = 0, and are located in the upper (lower) half-plane. 
The obvious generalization of Eqs. (Dll ) and (D12) in the case of r pairs of Bethe roots and holes is given by Eqs. (66) 
and (65 1. 



Appendix E: Derivation of the integral equations for the next-largest eigenvalues in the N/2 — 1 sector 

1. Integral equation for the auxiliary function 



As we have mentioned in Sec. IV C it is sufficient to consider the case with A^/2 — 1 Bethe roots and one hole in 



the relevant strip of the complexplane. A typical distribution of the Bethe roots and hole, at low-temperatures and 
T] e (0,7r/2), is presented in Fig. [3j where we have denoted by Xq the hole inside the strip |9fA| < rj/2. The eigenvalue 
and auxiliary function a s (A) corresponding to the distribution presented in Fig. [3] are described by formulas similar 
to (|47| and (48), but in this case q(X) is defined as 



9(A) = [] sinh(A - A,) , 



where {Xj}^l^ 1 are the A^/2 — 1 Bethe roots. The equation a s (A) + 1 = has 3V/2 — 1 solutions, of which, N/2 — 1 
are Bethe roots and V are holes denoted by {A^}^ 1 and A, 



Consider the contour C introduced in Sec. |IV A lj Inside the contour, the function 1 + a s (X) has N/2 zeros at the 
Bethe roots {A^}^^ 1 and hole Aq~, and a pole of order N/2 at iu. Therefore, we can define (A is located outside the 
contour C) 



fsW = tt~- / 7T (logsinh(A - /*)) log(l + a s (^))dfi = — / logsinh(A - n) , s[ dn , 
2m J c aX 2m J c l + a s (/x) 



which can be evaluated using Theorem [T] with the result 



N/2-1 



f s (X)= logsinh(A — Aj) + logsinh(A — Aq ) — — log sinh(A — iu) . 

j=i 



(El) 



(E2) 
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FIG. 3. Typical distribution of Bethe roots (•) and holes (o) in the strip I^Aj < 77 , 77 £ [0, 7r/2) characterizing the next-largest 
eigenvalues of the QTM in the N/2 — 1 sector without "particle-hole" type terms. All the other roots and holes can be obtained 
using the in periodicity. The contour C surrounds all the Bethe roots, the hole located at and the pole of the auxiliary 
function o s (A) at iu. 



Taking the logarithm in Eq. (48), and using (E2|, which can be rewritten as logg(A) = / S (A) — logsinh(A — A[j~) + 
log <j}-(^) + N/ 2 log sin 77 we find 



log a s ( A) = —(3h + log 



+ (A) 4>-{\ + iii) 
<t>-{\) <t>+(\ + ir)) 



log 



sinh(A — Aq — irj) 
sinh(A — Aq" + irj) 



+ f s (\ + ir,)-f s (\-ir,). 



Making use of Eq. (51), we can take the Trotter limit, N —> 00, obtaining the NLIE for the auxiliary function 

1 



log o fl (A) - -/3e (A + irj/2) t *tt + i0{\ - A+) 



2tt 



K(\-p) log(l + a s (/j))d/_t. 



(E3) 
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In Eq. ( E3 1 , the minus (plus) sign in front of the in factor is considered when Xq is in the first (second) quadrant of 
the complex plane. The same NLIE is valid also for r/ G (tt/2, n) if the contour C is replaced by a similar rectangular 
contour with the upper (lower) edges parallel to the real axis through ±i(n — r\ — e)/2 , e — > 0. 



2. Integral expression for the next-largest eigenvalue in the N/2 — 1 sector 



The starting point of our derivation will be again the representation (|55j), which is also valid for the eigenvalues in 
the N/2 — 1 sector. We need an integral representation for q^> (A) = sinh(A — Xq) Yif^ sinh(A — Ay ) . If we consider 
the contour C introduced in IV A 2[ then the following identity holds 



d(X-u) a ' s ^} du = Q. 
c+c> W l + a s (Ai) 



(E4) 



For A close to the real axis, using |E4|, Theorem [T] and the fact that inside the contour C' ph the function 1 + a s (A) 

has: N — 1 zeros at {A^ 
at — iu — if] + in, we find 



has: N - 1 zeros at {A^}^" 1 or { A^ + wr , N/2 - 1 simple poles at {Ay}^ -1 + if] and a pole of order N/2 



— J c d(X M ) l d^ = - f J2 d(A - Af) - J] dCA-A,--*,)--^ + »« + «?) 



(E5) 



Inside the contour C the function 1 + ft s (A) has N/2 zeros at the Bethe roots {AyVyVj 1 and hole Xq and a pole of 
order N/2 at iu. Using again Theorem [I] we have 



' ^ d(X — fi — ir})- ^ 



2ni 



N/2-1 



dji = d(X — Xj — irj) + d(X — Xq — iff) — ^rd{X — iu — it]) . 



1 + a s {u) 



(E6) 



Taking the difference of Eqs. (E5| and (E6|, integrating by parts, and then integrating w.r.t. A we obtain the following 
representation 



log gC*> (A) = log ( . , + A ° +) r ) + 1"M • ( A + , „ )o_ ( A / ,/ ) l 

,sinn(A — A — 177] , 



2tti 



[d(A - fi)-d{\- n - irj)} log(l + a s (u))d]i + c , 



(E7) 



with c a constant of integration. Th e integral expression for the next-largest eigenvalue of the QTM in the sector 
N/2 — 1 is obtained by replacing (E7) in (55) with the result 



, a Ph , ( sinhA+ A 1 

log A s = + log . + ° , + — . 

2 \smh(Ao + ir?)/ 2tt J c 



p {fi + ir]/2) log(l + a s [n))d(i. 



(E8) 



Eq. (E8| is also valid in the domain rj 6 (n/2,n) if the contour C is replaced by a rectangular contour, extending to 
infinity, with the edges parallel to the real axis through ±i(n — rj — e)/2 , e — > 0. 



3. Final form of the integral equations 



Performing the same operations as in Sec. IV A 3 and Appendix D3 Eq. (E3) is transformed into 



log o s (A - irj/2) = -Pe (X) T wr + i0(A - A+) 



1 

2n 



K(X - u) log(l + a s (u - ir}/2))dft , 



(E9) 



where we took into account that Aq — > Aq — i-q/2 The expression (E8) for the next-largest eigenvalue becomes 



logA s (0) 



/3/i 



*Po(Aj 



2tt 



(E10) 



3G 



Introducing the function v(X) satisfying e "( A )/ T = a s (X — irj/2) the NLIE for the auxiliary function and the integral 
expression for the next-largest eigenvalues in the N/2 — 1 sector at low-temperatures can be written as 



v(X) = e (A) ± inT - iT9(X - A+) + ^- [ K(X - p) log f 1 + e ~ v ^/ T \ dp, , 



(Ell) 



log A s (0) = — - in + ip Q (X^ 



1 

2n 



o(M)log( 



1 + e 



-v(m)/T 



dfi. 



(E12) 



where Xt satisfies the constraint 1 + e "( A o )/ T = 0. We should mention that we can discard the in term on the 



r.h.s. of ( E12 1 (this has the effect of neglecting an (— l) m factor in the asymptotic expansion which is irrelevant in 



the continuum limit) On the r.h.s. of Eq. (Ell) we will consider the plus (minus) sign in fro nt of the inT term when 
Xq is in the first (second) quadrant of the complex plane. The generalization of (Ell) and (E12) to the case when r 
"particle- hole" pairs are present is given by Eqs. (67) and (68). 



Appendix F: Proof of some identities 

Here we prove some identities used in Sec. [V] We start with 

p' (X)R(X,±q) dX=p' (q)-2np(q). 



(Fl) 



Using a formal solution of Eq. (33) on the l.h.s. of (Fl ) we have 



p' (X)R(X,±q)d\ 



q f ( 1 + ^ K K ) ( A < A ^ K (t* T 9)Po(A) d^ dX , 

Q 

-i 

where in the second line we have used the symmetry of the kernel K(X — /i) = K{[i — X) and the integral equation for 
the density (30). The identity (Fl) follows from 



1 f q 1 
p(±q) + 7^ J K(qT M)P(M) dn = ^Po( ±( l) 



and the fact that p(A) and p'o(X) are even functions. Using a similar method we can prove that 



e (X)p' (X)dX = 2n e (X)p(X)dX. 

q J-q 



(F2) 



(F3) 



Making use of the equation for the dressed energy (31) we can rewrite the l.h.s. of (F3) as 



2vr 



' q s (X)p' (X)dX= f p (l + ^-K 

' —q J —q J —q 

2n j e (p,)p(n)dfi, 

J-q 



(A,m) eoi^pQ^dpidX, 



where we have used again the symmetry of the kernel and Eq. ( 30 ) 



Appendix G: Derivation of the asymptotic expansion for the generating functional of density correlators 



In this Appendix we will show how we can derive using our method the results obtained in [54, 55 . The first step in 
the computation of the asymptotic expansion for the generating functional of density correlators in the Bose gas is the 
deri vation of NLIEs for the eigenvalues of the twisted QTM t®™ = A®™ (0) + e^D^™(0) in the N/2 sector (see 
Eq. (C21 ). The twisted QTM is a particular case of the generalized model considered in Appendix |B| with parameters 
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o(A) = b{u', \) N / 2 e f3h/2 , d(X) = b(X, - u ') N / 2 e -^ h / 2+v which means that the eigenvalues in the N/2 sector have the 
form 



a (¥>)/- \ a _ ^-(A) q(X-iy) jh 



0+(A) g(A + i?7)^_M 



^+(A + m7) e(A) 



with g(A) = Il^i smn (A — Aj) and {Aj}^f satisfying the BAEs 



b{u',\ 3 ) 



N/2 



_ e -l3h+ v 



N/2 

n 



sinh(A 3 — Afc + 177) 
sinh(Aj — A^. — 177) 



j = l,>--,N/2. 



(Gl) 



(G2) 



The derivation of the NLIEs and integral expressions for the N/2 sector eigenvalues of the twisted QTM (this includes 
also the largest eigenvalue) is almost identical with the one presented in Section IV A and Appendix [d] (we use a 
simil ar distribution of Bethe roots and holes as in F ig. | T| and Fig. [2]). We obtain similar equations as Eqs. (64a) and 
(64b) (for the largest eigenvalue) and Eqs. (65) and (|66| (for the next-largest eigenvalues in the N/2 sector) with the 
only difference being the replacement of the h]{2T) with h/(2T) + ipT in the r.h.s. of Eqs. (64b) and (66). The reader 
should note that ip does not appear in the integral expressions for the eigenvalues. 
Performing the continuum limit in (C21) we find 



00 . 



(G3) 



with the correlation lengths defined by 



f 



1 

'2tt 



log 



' 1 + e -«f(fc)/T > 



dk 



1 + e-^ k )/ T 

and the auxiliary functions uf(k) satisfying the NLIEs: 

uf{k) = k 2 - Li-<pT + iT^9(k - ifet) - iT^9(k -kj)-^ J K(k- k') log (l + e^^')/^ dk' . 



3=1 



3=1 



(G4) 



(G5) 



The 2r parameters, {k^Yj=i ({% }j=i) appearing in Eq. (G5) belong to the upper (lower) half of the complex plane 

and satisfy the constraint 1 + e Ui ^ T = 0. r can take the values 0, 1, 2, • • • with the r = ter m (this means that 
the sums in Eq. (G4), (G5) are zero) being the dominant contribution in the expansion. Eq. (G3) was first derived in 



Appendix H: Low-temperature limit of the asymptotic expansions 



Sec. 
and 



T he low - te mper ature analysis of the asymptotic expansions ( 10 1 and ( 13 1 is very similar with the one performed in 

1 — 



VB 



and 



vc 

k w 



The only difference is the fact that in the Bose gas case the principal integral operator is / — 



lich means that p' (k) = 1. Therefore, the calculations are almost identical with the ones for the 
XXZ spin chain except for some sign changes. The integral equations for the zero temperature dressed energy Eo(k) 
and the dressed charge Z(k) are given by 



e (k) 



1 f q ^ 



2tt 



K(k-k')e (k')dk' 



M = eo(fc) ■ 



(HI) 



and 



Z(k) 



2tt 



K(k - k')Z(k') dk' = 1 , Z(±q) = Z . 



(H2) 



The resolvent of the integral operator I — and the dressed phase equations are obtained from the XXZ spin 
chain equivalents, ( 33 ) and ( 34 ), by ch ang ing the sign in front of the integral and replacing K(\, [i) , 9(A) and ±g with 
K(k,k') ,9(k) and ±q. The identities (35) transform into (note the sign changes) 



Z(k) = 1 + F(k\ -q)- F(k\q) , = = 1 - F(q\q) - F(q\ - q) , 



(H3) 
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and ( Fl ) and ( F3 ) are still valid in the Bose case. Additional simplifications occur due to the fact that Z(k) = 2irp(k). 



The Fermi velocity can be rewritten as 



v F = 



4(<?) _ 



(H4) 



Using these relations and performing calculations similar with the ones from Sec. VB and VC we obtain (100) and 
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